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Abstract 

A measurement of inclusive W and Z production cross sections in pp collisions 
at v^ = 7TeV is presented. The electron and muon decay channels are ana- 
lyzed in a data sample collected with the CMS detector at the LHC and corre- 
sponding to an integrated luminosity of 36 pb^ . The measured inclusive cross 
sections are a{pp -^ WX) x 6(W -^ iv) = 10.30 ± 0.02 (stat.) ± 0.10 (syst.) ± 
0.10 (th.) ± 0.41 (lumi.) nb and cr{pp -^ ZX) x B{Z -^ £+£-) = 0.974 ± 
0.007 (stat.) ± 0.007 (syst.) ± 0.018 (th.) ± 0.039 (lumi.) nb, limited to the dilep- 
ton invariant mass range 60 to 120 GeV. The luminosity-independent cross sec- 
tion ratios are {cr{pp -^ WX) x B(W -^ iv)) / (c7(pp -^ ZX) x B{Z -^ £+£-)) = 
10.54 ±0.07 (stat.) ±0.08 (syst.) ±0.16 (th.) and {a{pp -^ W+X) x B(W+ -^ £+v)) / 
{a{pp -^ W-X) X B{W- -^ £-v)) = 1.421 ±0.006 (stat.) ±0.014 (syst.) ±0.029 (th.). 
The measured values agree with next-to-next-to-leading order QCD cross section cal- 
culations based on recent parton distribution functions. 
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1 Introduction 

This paper describes a measurement carried out by the Compact Muon Solenoid (CMS) Col- 
laboration of the inclusive production cross sections for W and Z bosons in pp collisions at 
^s = 7TeV. The vector bosons are observed via their decays to electrons and muons. In 
addition, selected cross-section ratios are presented. Precise determination of the production 
cross sections and their ratios provide an important test of the standard model (SM) of particle 
physics. 

The production of the electroweak (EWK) gauge bosons in pp collisions proceeds mainly via 
the weak Drell-Yan (DY) process lU consisting of the annihilation of a quark and an antiquark. 
The production process pp — ?■ W + X is dominated by ud — ?► W+ and du — ?► W^, while 
pp — )■ Z + X is dominated by uu and dd — ?► Z. 

Theoretical predictions of the total W and Z production cross sections are determined from 
parton-parton cross sections convolved with parton distribution functions (PDFs), incorpo- 
rating higher-order quantum chromodynamics (QCD) effects. PDF uncertainties, as well as 
higher-order QCD and EWK radiative corrections, limit the precision of current theoretical pre- 
dictions, which are available at next-to-leading order (NLO) [2 -4J and next-to-next-to-leading 
order (NNLO) USHU in perturbative QCD. 

The momentum fractions of the colliding partons Xj, X2 are related to the vector boson masses 
(^w/z ~ SX1X2) and rapidities {y = \ ln(xi/x2)). Within the accepted rapidity interval, \y\ < 
2.5, the values of x are in the range 10^^ < x < 0.1. 

Vector boson production in proton-proton collisions requires at least one sea quark, while two 
valence quarks are typical of pp collisions. Furthermore, given the high scale of the process, 
s = wj^/2 ~ 10^ GeV , the gluon is the dominant parton in the proton so that the scattering 
sea quarks are mainly generated by the g — > qq splitting process. For this reason, the precision 
of the cross section predictions at the Large Hadron Collider (LHC) depends crucially on the 
uncertainty in the momentum distribution of the gluon. Recent measurements from HERA IITOl 
and the Tevatron [.11.-.19J reduced the PDF uncertainties, leading to more precise cross-section 
predictions at the LHC. 

The W and Z production cross sections and their ratios w^ere previously measured by ATLAS [ [20| 
wth an integrated luminosity of 320 nb^^ and by CMS ||2T| with 2.9 pb^^. This paper presents 
an update with the full integrated luminosity recorded by CMS at the LHC in 2010, correspond- 
ing to 36 pb^^. The leptonic branching fraction and the width of the W boson can be extracted 
from the measured W/Z cross section ratio using the NNLO predictions for the total W and Z 
cross sections and the measured values of the Z boson total and leptonic partial widths ||22|. 
together wth the SM prediction for the leptonic partial width of the W 

This paper is organized as follows: in Section |2] the CMS detector is presented, with partic- 
ular attention to the subdetectors used to identify charged leptons and to infer the presence 
of neutrinos. Section |3] describes the data sample and simulation used in the analysis. The 
selection of the W and Z candidate events is discussed in Section Hj Section |5] describes the 
calculation of the geometrical and kinematic acceptances. The methods used to determine the 
reconstruction, selection, and trigger efficiencies of the leptons within the experimental accep- 
tance are presented in Section [6J The signal extraction methods for the W and Z channels, 
as well as the background contributions to the candidate samples, are discussed in Sections [7| 
and|8] Systematic uncertainties are discussed in Section|9] The calculation of the total cross sec- 
tions, along with the resulting values of the ratios and derived quantities, are summarized in 
SectionllO] In the same section w^e also report the cross sections as measured w^ithin the fiducial 
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and kinematic acceptance (after final-state QED radiation corrections), thereby eliminating the 
PDF uncertainties from the results. 



2 The CMS Detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, an electromagnetic calorimeter (ECAL), and a hadron calorimeter (HCAL). Muons are 
detected in gas-ionization detectors embedded in the steel return yoke. In addition to the barrel 
and endcap detectors, CMS has extensive forward calorimetry 

A right-handed coordinate system is used in CMS, with the origin at the nominal interaction 
point, the x-axis pointing to the center of the LHC ring, the y-axis pointing up (perpendicular 
to the LHC plane), and the z-axis along the anticlockwise-beam direction. The polar angle 6 
is measured from the positive z-axis and the azimuthal angle cp is measured (in radians) in the 
xy-p\ane. The pseudorapidity is given by ^ = — lntan(0/2). 

The inner tracker measures charged particle trajectories in the pseudorapidity range \^\ < 2.5. 
It consists of 1440 silicon pixel and 15 148 silicon strip detector modules. It provides an impact 
parameter resolution of ~15 }im and a transverse momentum (pj) resolution of about 1% for 
charged particles with pj ^ 40 GeV. 

The electromagnetic calorimeter consists of nearly 76 000 lead tungstate crystals, which provide 
coverage in pseudorapidity |?/| < 1.479 in a cylindrical barrel region (EB) and 1.479 < \i]\ < 3.0 
in two endcap regions (EE). A preshower detector consisting of two planes of silicon sensors 
interleaved with a total of three radiation lengths of lead is located in front of the EE. The ECAL 
has an energy resolution of better than 0.5% for unconverted photons with transverse energies 
(Ej) above 100 GeV. The energy resolution is 3% or better for the range of electron energies rel- 
evant for this analysis. The hadronic barrel and endcap calorimeters are sampling devices with 
brass as the passive material and scintillator as the active material. The combined calorimeter 
cells are grouped in projective towers of granularity Arj x A^ = 0.087 x 0.087 at central rapidi- 
ties and 0.175 x 0.175 at forward rapidities. The energy of charged pions and other quasi-stable 
hadrons can be measured w^ith the calorimeters (ECAL and HCAL combined) with a resolution 
of AE/E ~ 100%/y E( GeV) © 5%. For charged hadrons, the calorimeter resolution improves 
on the tracker momentum resolution only for pj in excess of 500 GeV. The energy resolution 
on jets and missing transverse energy is substantially improved with respect to calorimetric 
reconstruction by using the particle flow^ (PF) algorithm fZ\ which consists in reconstructing 
and identifying each single particle with an optimised combination of all sub-detector infor- 
mation. This approach exploits the very good tracker momentum resolution to improve the 
energy measurement of charged hadrons. 

Muons are detected in the pseudorapidity window \r]\ < 2.4, with detection planes based on 
three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. A high-px 
muon originating from the interaction point produces track segments typically in three or four 
muon stations. Matching these segments to tracks measured in the inner tracker results in a px 
resolution between 1 and 2% for px values up to 100 GeV. 



The first level (LI) of the CMS trigger system ||24J, composed of custom hardware processors, 
is designed to select the most interesting events in less than 1 fis, using information from the 
calorimeters and muon detectors. The High Level Trigger (HLT) processor farm p5| further 
decreases the event rate to a few hundred Hz before data storage. A more detailed description 



of CMS can be found elsewhere 

3 Data and Simulated Samples 

The W and Z analyses are based on data samples collected during the LHC data operation 
periods logged from May through November 2011, corresponding to an integrated luminosity 
/lint = 35.9 ± 1.4 pb"\ 

Candidate events are selected from datasets collected with high-Ej lepton trigger requirements. 
Events w^ith high-Ex electrons are selected online if they pass a LI trigger filter that requires an 
energy deposit in a coarse-granularity region of the ECAL with £x > 5 or 8 GeV, depending on 
the data taking period. They subsequently must pass an HLT filter that requires a minimum Ej 
threshold of the ECAL cluster which is well below the offline Ej threshold of 25 GeV. The full 
ECAL granularity and offline calibration corrections are exploited by the HLT filter |j27|. 

Events with high-px muons are selected online by a single-muon trigger. The energy threshold 
at the LI is 7 GeV. The px threshold at the HLT level depends on the data taking period and 
was 9 GeV for the first 7.5 pb^^ of collected data and 15 GeV for the remaining 28.4 pb^^. 

Several large Monte Carlo (MC) simulated samples are used to evaluate signal and background 
efficiencies and to validate the analysis techniques employed. Samples of EWK processes with 
Z and W bosons, both for signal and background events, are generated using POWHEG II28I - 
|30l interfaced with the PYTHIA |31| parton-shower generator and the Z2 tune (the PYTHIA6 
Z2 tune is identical to the Zl tune described in I.32J except that Z2 uses the CTEQ6L PDF, 
while Zl uses the CTEQ5L PDF). QCD multijet events with a muon or electron in the final 
state and tt events are simulated with pythia. Generated events are processed through the full 
Geant4 Il33ll34| detector simulation, trigger emulation, and event reconstruction chain of the 
CMS experiment. 

4 Event Selection 

The W — > £t/ events are characterized by a prompt, energetic, and isolated lepton and signifi- 
cant missing transverse energy, ^j. No requirement on ^j is applied. Rather, the ^x is used as 
the main discriminant variable against backgrounds from QCD events. 

The Z boson decays to leptons (electrons or muons) are selected based on two energetic and 
isolated leptons. The reconstructed dilepton invariant mass is required to be consistent with 
the known Z boson mass. 

The following background processes are considered: 

• QCD multijet events. Isolation requirements reduce events with leptons produced 
inside jets. The remaining background is estimated with a variety of techniques 
based on data. 

• High-Ej photons. For the W — > et/ channel only, there is a nonnegligible back- 
ground contribution coming from the conversion of a photon from the process pp — >■ 
7+jet(s). 

• Drell-Yan. A DY lepton pair constitutes a background for the W — ?► £i/ channels 
when one of the two leptons is not reconstructed or does not enter a fiducial region. 

• W — 7- Ti/ and Z — ?► T+T^ production. A small background contribution comes from 
W and Z events with one or both t decaying leptonically The minimum lepton px 
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requirement tends to suppress these backgrounds. 

• Diboson production. The production of boson pairs (WW, WZ, ZZ) is considered 
a background to the W and Z analysis because the theoretical predictions for the 
vector boson production cross sections used for comparison with data do not include 
diboson production. The background from diboson production is very small and is 
estimated using simulations. 

• Top-quark pairs. The background from tt production is quite small and is estimated 
from simulations. 

The backgrounds mentioned in the first two bullets are referred to as "QCD backgrounds", the 
Drell-Yan, W — > rv, and dibosons as "EWK backgrounds", and the last one as "tt background". 
For both diboson and tt backgrounds, the NLO cross sections were used. The complete selec- 
tion criteria used to reduce the above backgrounds are described below. 

4.1 Lepton Isolation 

The isolation variables for the tracker and the electromagnetic and hadronic calorimeters are 
defined: Jtrk = Etracks Pt , kcAL = Eecal ^t , Jhcal = Ehc al ^t , where th e sums are per- 
formed on all objects falling w^ithin a cone of aperture AR = y^(Af/)2 + (A^)^ = 0.3 around the 
lepton candidate momentum direction. The energy deposits and the track associated wth the 
lepton candidate are excluded from the sums. 

4.2 Electron Channel Selection 

Electrons are identified offline as clusters of ECAL energy deposits matched to tracks recon- 
structed in the silicon tracker. The ECAL clustering algorithm is designed to reconstruct clus- 
ters containing a large fraction of the energy of the original electron, including energy radiated 
along its trajectory. The ECAL clusters must fall in the ECAL fiducial volume of |?/| < 1.44 
for EB clusters or 1.57 < |?/| < 2.5 for EE clusters. The transition region 1.44 < \rj\ < 1.57 is 
excluded as it leads to lower-quality reconstructed clusters, due mainly to services and cables 
exiting between the barrel and endcap calorimeters. Electron tracks are reconstructed using an 
algorithm [35 1 (Gaussian-sum filter, or GSF tracking) that accounts for possible energy loss due 
to bremsstrahlung in the tracker layers. 

The radiated photons may convert close to the original electron trajectory, leading to charge 
misidentification. Three different methods are used to determine the electron charge. First, the 
electron charge is determined by the signed curvature of the associated GSF track. Second, the 
charge is determined from the associated trajectory reconstructed in the silicon tracker using a 
Kalman Filter algorithm Il36l . Third, the electron charge is determined based on the azimuthal 
angle betw^een the vector joining the nominal interaction point and the ECAL cluster position 
and the vector joining the nominal interaction point and innermost hit of the GSF track. The 
electron charge is determined from the two out of three charge estimates that are in agreement. 
The electron charge misidentification rate is measured in data using the Z — ?> e+e^ data sample 
to be within 0.1%-1.3% in EB and 1.4%-2.1% in EE, increasing w^ith electron pseudorapidity 

Events are selected if they contain one or tw^o electrons having Ej > 25 GeV for the W — > ev 
or the Z — )■ e+e^ analysis, respectively. For the Z — > e+e^ selection there is no requirement 
on the charges of the electrons. The energy of an electron candidate with Ej > 25 GeV is 
determined by the ECAL cluster energy, while its momentum direction is determined by that 
of the associated track. 

Particles misidentified as electrons are suppressed by requiring that the rj and (p coordinates of 
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Figure 1: Distributions of the electron identification variables At], Acp, cr,^,^, and H/E for data 
(points with the error bars), for EB (left) and EE (right). For illustration the simulated W — > ev 
signal (histograms), normalized to the number of events observed in data, is superimposed. 
These distributions are obtained after applying all the tight requirements on the selection vari- 
ables, except that on the presented variable. The tight requirement on that variable is indicated 
with an arrow. 
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Figure 2: Distributions of the electron isolation variables Itj-k/Ei/ Jecal/^T/ arid /hcal/^t 
for data (points with the error bars), for EB (left) and EE (right). For illustration the simu- 
lated W — > et/ signal (histograms), normalized to the number of events observed in data, is 
superimposed. These distributions are obtained after applying all the tight requirements on 
the selection variables, except that on the presented variable. The tight requirement on that 
variable is indicated with an arrow. 
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the track trajectory extrapolated to the ECAL match those of the ECAL cluster permitting only 
small differences {Atj, A(p) between the coordinates, by requiring a narrow ECAL cluster width 
in ?! {cTfjrj), and by limiting the ratio of the hadronic energy H to the electromagnetic energy 
E measured in a cone of AR = 0.15 around the ECAL cluster direction. More details on the 
electron identification variables can be found in Refs. ||37ll38l| . Electron isolation is based on 
requirements on the three isolation variables Ihcal/Et/ J^ecal/Et/ arid /trk/Ej- 

Electrons from photon conversions are suppressed by requiring the reconstructed electron track 
to have at least one hit in the innermost pixel layer. Furthermore, electrons are rejected when a 
partner track is found that is consistent with a photon conversion, based on the opening angle 
and the separation in the transverse plane at the point w^here the electron and partner tracks 
are parallel. 

The electron selection criteria were obtained by optimizing signal and background levels ac- 
cording to simulation-based studies. The optimization was done for EB and EE separately. 

Two sets of electron selection criteria are considered: a tight one and a loose one. Their effi- 
ciencies, from simulation studies based on W — > ei/ events, are approximately 80% and 95%, 
respectively. These efficiencies correspond to reconstructed electrons w^ithin the geometrical 
and kinematic acceptance, which is defined in Section |5] The tight selection criteria give a 
purer sample of prompt electrons and are used for both the W — > ev and Z — ?> e+e^ analyses. 
The virtue of this choice is to have consistent electron definitions for both analyses, simplifying 
the treatment of systematic uncertainties in the W/Z ratio measurement. In addition, the tight 
working point, applied to both electrons in the Z — > e+e^ analysis, reduces the QCD back- 
grounds to a negligible level. Distributions of the selection variables are shown in Figs.[l]and|2] 
The plots show the distribution of data together with the simulated signal normalized to the 
same number of events as the data, after applying all the tight requirements on the selection 
variables except the requirement on the displayed variable. 

For the W analysis, an event is also rejected if there is a second electron that passes the loose 
selection with Ej > 20 GeV. This requirement reduces the contamination from DY events. 
The number of W — ?► ev candidate events selected in the data sample is 235 687, with 132 696 
positrons and 102 991 electrons. 

For the Z analysis, two electrons are required within the ECAL acceptance, both with Ej > 
25 GeV and both satisfying the tight electron selection. Events in the dielectron mass region of 
60 < ffJee < 120 GeV are counted. These requirements select 8452 events. 

4.3 Muon Channel Selection 

Muons candidates are first reconstructed separately in the central tracker (referred to simply 
as "tracks" or "tracker tracks") and in the muon detector ("stand-alone muons"). Stand-alone 
muons are then matched and combined with tracker tracks to form "global muons". Another 
independent algorithm proceeds from the central tracker outwards, matching muon chambers 
hits and producing "tracker muons". 

The following quality selection are applied to muon candidates. Global and stand-alone muon 
candidates must have at least one good hit in the muon chambers. Tracker muons must match 
to hits in at least two muon stations. Tracks, global muons, and tracker muons must have 
more than 10 hits in the inner tracker, of w^hich at least one must be in the pixel detector, and 
the impact parameter in the transverse plane, dxy, calculated with respect to the beam axis, 
must be smaller than 2 mm. More details and studies on muon identification can be found in 
Refs. 1391 iQl. 
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Figure 3: Distribution of number of hits in the inner tracker and in the pixel detector, num- 
ber of hits in muon chambers, number of muon segments stations, x^ per degree of freedom, 
and transverse impact parameter dxij for data (points with the error bars). For illustration the 
simulated W — ?► /^v signal (histogram), normalized to the number of events observed in data, 
is superimposed. These distributions are for events satisfying all selection requirements, ex- 
cept that on the presented variable. The applied requirement on that variable is indicated as 
a blue arrow. In the dxy distribution, the horizontal line shows the average of the bins with 
dxy > 0.2 cm used to estimate the cosmic-ray muon contamination in the signal region. The ex- 
cess of events in data in the region with d^y < 0.2 cm with respect to W — ?► f/i/ signal simulation 
is due to muons from long-lived particle decays in the QCD background. 



4.3 Muon Channel Selection 



Muon candidates selected in the W — ?► /^v analysis must be identified both as global and tracker 
muons. Moreover, as additional quality selection, the global muon fit must have a }^ per degree 
of freedom less than 10 in order to reject misidentified muons and misreconstructed particles. 
The W — ?■ /^i/ candidate events must have a muon candidate in the fiducial volume |//| < 2.1 
with Pt > 25 GeV. The muon must be isolated, satisfying JJ^J^^^ = (Jt^k + ^ecal + ^hcal) /Pt < 
0.1. Events containing a second muon with px > 10 GeV in the full muon acceptance region 
(I?/ 1 < 2.4) are rejected to minimize the contamination from DY events. The distributions of 
the variables used for muon quality selection are shown in Fig. [3] after applying all selection 
requirements, except that on the presented variable. 

Background due to a cosmic-ray muon crossing the detector in coincidence with a pp colli- 
sion is very much reduced by the impact parameter requirement. The remaining cosmic-ray 
background is evaluated by extrapolating to the signal region the rate of events with large im- 
pact parameter. Figure [3| (bottom, right) shows the distribution of the impact parameter dxy for 



the W — > /^i/ candidates satisfying all selection requirements, except that on d^y. Candidates 
with large d^y are mainly due to cosmic-ray muons and their rate is independent of dxy- A 
background fraction on the order of 10^^ in the dxy < 2 mm region is estimated. 

The isolation distribution in data, together with the MC expectations, are shown in Fig. HI 
Events with I'^f^^ > 0.2 are mainly from QCD multijet background, and are used as a control 
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Figure 4: Distribution of Ilf^^ for candidates with a good quality muon of pj > 25 GeV in the 
fiducial region \t]\ < 2.1. Points represent the data and the histograms the contribution from 
the different SM processes. The signal selection requirement (dark blue arrow, Ilf^^ < 0.1) and 
the selection of the QCD-enriched control sample (light green arrowy, Ilf^^ > 0.2) are shown. 

sample (Section [73) |. 

After the selection process described, 166457 events are selected, 97533 of them with a posi- 
tively charged muon candidate and 68 924 with a negatively charged muon candidate. 

Z — )■ n^n^ candidate events are selected by pairing a global muon matched to an HLT trigger 
muon w^ith a second oppositely charged muon candidate that can be either a global muon, a 
stand-alone muon, or a track. No ]^ selection or requirement that the muon be reconstructed 
through the tracker-muon algorithm is applied. The tw^o muon candidates must both have 
px > 20 GeV and \r]\ < 2.1, and their invariant mass must be in the range 60 < m^^ < 120 GeV. 
Both muon candidates must be isolated according to the tracker isolation requirement 7trk < 
3 GeV. The different choice of isolation requirements in W — ?► }iv and Z — ?► f/+/i^ is motivated in 



Section 8.3 After the selection process, the number of selected events with two global muons 
is 13 728. 



10 



6 Efficiencies 



5 Acceptance 

The acceptance Aw(e) for W — )■ ev is defined as the fraction of simulated W events having an 
ECAL cluster within the ECAL fiducial volume with Ej > 25 GeV. The ECAL cluster must 
match the generated electron after final-state radiation (FSR) within a cone of AR = 0.2. No 
matching in energy is required. 

There is an inefficiency in the ECAL cluster reconstruction for electrons direction within the 
ECAL fiducial volume due to a small fraction (0.5%) of noisy or malfunctioning towers re- 
moved from the reconstruction. These are taken into account in the MC simulation, and no 
uncertainty is assigned to this purely geometrical inefficiency. The ECAL cluster selection effi- 
ciency is also affected by a bias in the electron energy scale due to the 25 GeV energy threshold. 
The related systematic uncertainty is assigned to the final W and Z selection efficiencies. 

The acceptance for the Z — > e+e^ selection, Az{e), is defined as the number of simulated 
events with tw^o ECAL clusters with Ej > 25 GeV w^ithin the ECAL fiducial volume and with 
invariant mass in the range 60 < nigg < 120 GeV, divided by the total number of signal events 
in the same mass range, with the invariant mass evaluated using the momenta at generator 
level before FSR. The ECAL clusters must match the two simulated electrons after FSR within 
cones of AR < 0.2. No requirement on energy matching is applied. 

For the W ^ }iv analysis, the acceptance Aw (fi) is defined as the fraction of simulated W signal 
events w^ith muons having transverse momentum pj^'^ and pseudorapidity ?/§^", evaluated at 
the generator level after FSR, within the kinematic selection: pl*^" > 25 GeV and If/S*^"! < 2.1. 

The acceptance Az{}i) for the Z — ?> fi^ji^ analysis is defined as the number of simulated Z sig- 
nal events with both muons passing the kinematic selection with momenta evaluated after FSR, 
Pj*^" > 20 GeV and |f/8en| ^ 2.I, and with invariant mass in the range 60 < m^^ < 120 GeV, 
divided by the total number of signal events in the same mass range, with the invariant mass 
evaluated using the momenta at generator level before FSR. 

Table n^ presents the acceptances for W+, W^, and inclusive W and Z events, computed from 
samples simulated w^ith POWHEG using the CTIO PDF, for the muon and the electron channels. 
The acceptances are affected by several theoretical uncertainties, which are discussed in detail 
in Section [931 

Table 1: Acceptances from POWHEG (with CTIO PDF) for W -> iv and Z -> £+£- final states, 
with the MC statistics uncertainties. 



Process 



Aw,2 



i = ^ 



W+ -^ i+v 0.5017 ± 0.0004 
W- -^ £-v 0.4808 ± 0.0004 

w ^ ev 0.4933 ± 0.0003 



0.4594 ± 0.0004 
0.4471 ± 0.0004 
0.4543 ± 0.0003 



i+t 



0.3876 ± 0.0005 



0.3978 ± 0.0005 



6 Efficiencies 

A key component of this analysis is the estimation of lepton efficiencies. The efficiency is de- 
termined for different selection steps: 
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• offline reconstruction of the lepton; 

• lepton selection, with identification and isolation criteria; 

• trigger (Ll+HLT). 

The order of the above selections steps is important. Lepton efficiency for each selection is 
determined with respect to the prior step. 

A tag-and-probe (T&P) technique is used, as described below, on pure samples of Z — > £+^^ 
events. The statistical uncertainty on the efficiencies is ultimately propagated as a systematic 
uncertainty on the cross-section measurements. This procedure has the advantage of extracting 
the efficiencies from a sample of leptons kinematically very similar to those used in the W 
analysis and exploits the relatively pure selection of Z — > £+^^ events obtained after a dilepton 
invariant mass requirement around the Z mass. 

The T&P method is as follows: one lepton candidate, called the "tag", satisfies trigger criteria, 
tight identification and isolation requirements. The other lepton candidate, called the "probe", 
is required to pass specific criteria that depend on the efficiency under study. 

For each kind of efficiency, the T&P method is applied to real data and to simulated samples, 
and the ratio of efficiencies in data (edata) arid simulation (egim) is computed: 

fcsim 

together with the associated statistical and systematic uncertainties. 

6.1 Electrons 

As mentioned in the previous section, the tight electron selection is considered for both the W 
and Z analyses, so: 

^all = ^rec^tight^trg- (2) 

The reconstruction efficiency ^rec is relative to ECAL clusters within the ECAL acceptance, the 
selection efficiency ^tight is relative to GSF electrons within the acceptance, and the trigger effi- 
ciency etrg is relative to electrons satisfying the tight selection criteria. 

All the efficiencies are determined by the T&P technique. Selections with different criteria have 
been tried on the tag electron. It was found that the estimated efficiencies are insensitive to 
the tag selection definition. The invariant mass of the T&P pair is required to be within the 
window 60 < m^g < 120 GeV, ensuring high purity of the probe sample. No opposite-charge 
requirement is enforced. 

The number of probes passing and failing the selection is determined from fits to the invariant 
mass distribution, with signal and background components. Estimated backgrounds, mostly 
from QCD multijet processes, are in most cases at the percent level of the overall sample, but 
can be larger in subsamples where the probe fails a selection, hence the importance of back- 
ground modeling. The signal shape is a Breit-Wigner w^ith nominal Z mass and width con- 
volved with an asymmetric resolution function (Crystal Ball HT)) with floating parameters. 
The background is modeled by an exponential. Systematic uncertainties that depend on the ef- 
ficiency under study are determined by considering alternative signal and background shape 
models. Details can be found in Section HI 

The T&P event selection efficiencies in the simulation are determined from large samples of 
signal events with no background added. 
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The T&P efficiencies are measured for the EB and EE electrons separately. Tag-and-probe effi- 
ciencies are also determined separately by charge, to be used in the measurements of the W+ 
and W^ cross sections and their ratio. Inclusive efficiencies and correction factors are summa- 
rized in Table |2] The T&P measurements of the efficiencies on the right-hand side of Eq. (|2]| are 
denoted as et&p-rec/ ^t&p-tight/ and et&p-trg- 

Table 2: Tag-and-probe efficiencies in data and simulation, and the correction factors used 
in the electron channels for the barrel (EB) and endcaps (EE). The combined statistical and 
systematic uncertainties are quoted. 



Efficiency Data 


Simulation 


Data /simulation (p) 


EB 


£t&p-rec 

^t&p-tight 

^t&p-trg 


(97.0 ± 1.0)% 

(84.0 ±0.3)% 
(98.0 ±0.1)% 


(97.78 ±0.02)% 
(87.47 ±0.05)% 
(97.10 ±0.03)% 


0.992 ±0.011 
0.960 ± 0.004 
1.009 ±0.001 


^t&p-aU 


(79.8 ±0.9)% 


(83.05 ± 0.06)% 


0.961 ±0.011 


EE 


^t&p-rec 

^t&p-tight 

^t&p-trg 


(94.3 ±1.1)% 
(73.1 ±0.7)% 
(97.3 ±0.3)% 


(94.61 ±0.05)% 
(75.61 ±0.06)% 
(97.16 ±0.04)% 


0.997 ±0.011 
0.966 ± 0.009 
1.001 ±0.003 


^t&p-aU 


(67.0 ± 1.0)% 


(69.51 ±0.07)% 


0.965 ± 0.015 



Event selection efficiencies are measured with respect to the W events within the ECAL ac- 
ceptance. Simulation efficiencies estimated from POWHEG W samples are shown in Table |3] 
These are efficiencies at the event level, e.g.: they include efficiency loss due to the second elec- 
tron veto. Given the acceptances listed in Table [l] and the T&P efficiencies listed in Table |2l 
the overall efficiency correction factors for electrons from W decays are computed. The over- 
all W signal efficiencies, obtained as products of simulation efficiencies with data/ simulation 
correction factors, are listed in Table |3] 

Table 3: Simulation efficiencies and the final corrected selection efficiencies for the W+, W^, 
and their average, in the W — > ev analysis. The quoted uncertainties are statistical for egim and 
include both statistical and systematic uncertainties for the corrected efficiencies egim x P- 





£sim 


^sim X p 


W+ -^ e+v 
W ^ ei/ 


(76.04 ±0.03)% 
(76.94 ±0.03)% 
(76.40 ±0.02)% 


(73.7 ±1.0)% 
(73.2 ± 1.0)% 
(73.5 ±0.9)% 



The efficiencies and the data/ simulation ratios are also estimated in bins of the electron Ex and 
rj in order to examine in detail the detector performance and take into account the differences in 
the W and Z kinematic distributions. The data /simulation ratios for reconstruction, selection, 
and trigger are shown in Fig. [5] as functions of the electron Ej and ?/. 

The reconstruction data /simulation ratios appear to be uniform with respect to Ej and rj, so a 
smaller number of bins is sufficient for the determination of their values. The data/simulation 
ratios for the selection and trigger efficiencies show a dependence that is estimated using ten ?/ 
bins and six Ej bins. Data /simulation ratios are estimated for both electron charges as well. 

The binned ratios and simulation efficiencies are transferred into the W analysis by properly 
weighting their product in each (Ej, f]) bin by the relative ECAL cluster abundance estimated 
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from POWHEG simulations. The corrected efficiencies are compared with the two-bin case in 
which the efficiencies are estimated in two bins of rj (EB and EE). The multibin corrected effi- 
ciencies are found to be consistent with the two-bin corrected efficiencies within the assigned 
uncertainties. In order to be sure that no hidden systematic uncertainty is missed, half of the 
maximum difference between the multibin and two-bin corrected efficiencies is propagated 
as an additional systematic uncertainty on the two-bin efficiencies used to estimate the cross 
sections. The additional relative uncertainty is at the level of 0.6%. 

The Z selection efficiencies for data and simulation are obtained based on the T&P efficiencies 
listed in Table [2] and the event acceptances given in Table IT] The Z efficiencies are first deter- 
mined after reconstruction and identification (as products of single-electron efficiencies). The 
event trigger efficiency is computed as the probability that at least one of the two electrons sat- 
isfies the Ll+HLT requirement. The overall selection efficiency for the Z analysis is the product 
of the reconstruction, identification, and trigger efficiencies. The simulation efficiency obtained 
from the POWHEG Z samples, together with the final corrected Z selection efficiency egim x p, 
are shown in Table H) These efficiencies are relative to the Z events with both electrons within 
the ECAL acceptance. 

Table 4: Simulation efficiency and the final corrected selection efficiency for the Z — ?> e+e^ 
analysis. The quoted uncertainties are statistical for egim and include both statistical and sys- 
tematic uncertainties for the corrected efficiency ^sim x P- 





Muons 




£'sim 


^sim X p 




Z -^ e+e- 


(66.74 ±0.07)% 


(60.9 ±1.1)% 


6.2 









For the W — )■ /^i/ cross section determination the single-muon efficiency combines the efficien- 
cies of all the steps in the muon selection: triggering on the muon, reconstructing it in the muon 
and central detectors, and applying the quality selection and the isolation requirement. In the 
procedure followed in this analysis, the reconstruction efficiency in the central tracker is factor- 
ized and computed independently, while the remaining terms are computed globally, without 
further factorizing them into different terms. 

An initial preselection of Z events for the T&P method is performed by selecting events that 
contain tracks measured in the central tracker having pj > 25 GeV, \f]\ < 2.1, and, when com- 
bined w^ith an oppositely charged track, give an invariant mass in the range 60 < m„+„- < 
120 GeV. We further require the presence in the event of a "tag" muon, defined as a global 
muon, that is matched to one of the preselected tracks, passes the selection described in Sec- 
tion 4.3 and corresponds to an HLT muon. The number of tag muons selected in data is about 
22 000. All the other preselected tracks are considered as probes to evaluate the muon efficiency. 
The background present in this sample is subtracted with a fit to the dimuon invariant mass 
spectrum of the sum of a Z component and a linear background contribution. The shape of the 
Z component is taken from simulation. 

The efficiency is studied as a function of the muon rj and pj. A dependence on rj is observed 
(Fig.|6J left) because different regions are covered by different muon detectors. This behavior is 
not fully reproduced in the simulation, as reflected in the corresponding p values (Fig. 16} right). 
The efficiency also exhibits a dependence on pj (Fig. U\ left), but this trend is similar in data 
and in simulation, and the correction factors can be taken as approximately constant up to px 
= 100 GeV (Fig. It] right). These binned correction factors are applied to the W analysis during 
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Figure 5: Data /simulation T&P ratios versus electron Ej (left column) and rj (right column). 
The ratios are presented for the reconstruction {prec, top row), selection (ptight/ middle row), 
and trigger (ptrg/ bottom row) efficiencies. Points with error bars represent the ratio measured 
in data; dashed lines correspond to a constant ratio of one. 
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Figure 6: Single-muon efficiencies (left) for data (red circles with error bars) ar\d simulation 
(black triangles), and the ratio between them (right), as a function of the muon tj. 
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Figure 7: Single-muon efficiencies (left) for data (red circles with error bars) and simulation 
(black triangles) and the ratio betw^een them (right), as a function of the muon pj. 

signal modeling (Section 8): W simulated events are w^eighted with the p factor corresponding 
to the (pt/ ?/) bin of the muon. The slightly difference between the kinematic characteristics of 
the muons and those from W decays is thus taken into account. 

The average single-muon efficiencies and correction factors are reported in Table |5] for posi- 
tively and negatively charged muons separately, and inclusively. The statistical uncertainties 
reflect the size of the available Z sample. Systematic uncertainties on et&p(data) and the cor- 



rection factors p are discussed in Section 9.2 



A small fraction of muon events are lost because of LI muon trigger prefiring, i.e., the assign- 
ment of a muon segment to an incorrect bunch crossing, occurring with a probability of a few 
per mille per segment. The effect is only sizable in the drift-tube system. The efficiency correc- 



16 



7 The W -> £i/ Signal Extraction 



Table 5: Tag-and-probe efficiencies in data and simulation and correction factors for positively 
and negatively charged muons. The errors on et&p (sim) are statistical only, while the systematic 
uncertainty is included for the other quantities. 



et&p (data) 
et&p (sim) 

P 



H 



¥ 



r 



(86.0 ±0.8)% 

(89.25 ±0.05)% 

(96.3 ±0.9)% 



(85.0 ±0.8)% 

(89.38 ±0.05)% 

(95.1 ±0.9)% 



(85.58 ±0.8)% 

(89.32 ±0.04)% 

(95.7 ±0.9)% 



tion in the barrel region is estimated for the current data to be ~1% per muon. This estimate is 
obtained from studies of muon pairs selected by online and offline single-muon trigger paths 
at the wrong bunch crossing, that have an invariant mass near the Z mass. Tracker information 
is not present in the case of prefiring, precluding the building of a trigger muon online or a 
global muon in the offline reconstruction. Since this effect is not accounted for in the efficiency 
from T&P, the measured Z — > ji^ji^ and W — > /^t/ cross sections are increased by 1% and 0.5%, 
respectively (including barrel and endcap regions) to correct for the effect of trigger prefiring. 
The uncertainty on those corrections is taken as a systematic uncertainty. 



The W — > /^i/ efficiencies from simulation are shown in Table ^ for the W+ and W^ samples 
separately and combined after applying the binned corrections estimated w^ith the T&P method 
using Z events. 

Table 6: Simulation efficiencies and final corrected efficiencies for the W — > /^v analysis. The 
quoted uncertainties are statistical for esim and include both statistical and systematic uncer- 
tainties for the corrected efficiencies egim x P- 




(89.19 ±0.03)% 
(89.19 ±0.03)% 
(89.19 ±0.03)% 



£sim X P 

(85.4 ±0.8)% 
(84.1 ±0.8)% 
(84.8 ±0.8)% 



For the Z ^ ji^fi cross section measurement, the muon efficiencies are determined together 



with the Z yield using a simultaneous fit described in Section 8.3 



7 The W ^>- £v Signal Extraction 

The signal and background yields are obtained by fitting the ^t distributions for W — )■ ev 
and W — 7- }iv to different functional models. An accurate I^j measurement is essential for 
distinguishing a W signal from QCD multijet backgrounds. We profit from the application of 
the PF algorithm, which provides superior ^j reconstruction performance ||42| w^ith respect to 
alternative algorithms at the energy scale of the W boson. 

The Sj is the magnitude of the transverse component of the missing momentum vector, com- 
puted as the negative of the vector sum of all reconstructed transverse momenta of particles 
identified with the PF algorithm. The algorithm combines the information from the inner 
tracker, the muon chambers, and the calorimeters to classify reconstructed objects according 
to particle type (electron, muon, photon, or charged or neutral hadron), thereby allowing pre- 
cise energy corrections. The use of the tracker information reduces the sensitivity of Sj to 
miscalibration of the calorimetry 

The QCD multijet background is one of the most significant backgrounds in W analyses. At 



7.1 Signal ^T Modeling 17 

high ^T/ EWK backgrounds, in particular W — )■ ti/ and DY, also become relevant, leading to 
contamination levels on the order of 10%. 

The Bt model is fitted to the observed distribution as the sum of three contributions: the W 
signal, and the QCD and EWK backgrounds. The EWK contributions are normalized to the W 
signal yield in the fit through the ratios of the theoretical cross sections. 

Simultaneous fits are performed to the two ^t spectra of W+ and W^ candidates, fitting ei- 
ther the total W cross section and the ratio of positive and negative W cross sections, or the 
individual positive and negative W cross sections. In both cases the overall normalization of 
QCD multijet events is determined from the fit. The diboson and tt contributions, taken from 
simulations, are negligible (Section |7.2| . 

In the following sections the modeling of the ^j shape for the signal and the EWK backgrounds 
are presented, and the methods used to determine the ^j shape for the QCD multijet back- 
ground from data are described. Finally, the extraction of the signal yields is discussed. 

7.1 Signal £r Modeling 

The W — ?> ^1/ signal is extracted with methods that employ simulation predictions of the Sj 
distribution in signal events. These predictions rely on the modeling of the vector-boson recoil 
and detector effects that can be difficult to simulate accurately. Discrepancies could result from 
deficiencies in the modeling of the calorimeter response and resolution, and from an incomplete 
description of the underlying event. These residual effects are addressed using corrections 
determined from the study of Z-boson recoil in data, discussed in the following paragraph. 

The recoil to the vector boson is defined as the negative of the vector sum of transverse energy 
vectors of all particles reconstructed with the PF algorithm in W and Z events, after subtract- 
ing the contribution from the daughter lepton(s). The recoil is determined for each event in 
Z — > £^£~ data and simulated Z — >• £~^£^ and W — > £v samples. We fit the distributions of 
the recoil components (parallel and perpendicular to the boson p^ direction) with a double 
Gaussian, whose mean and width vary with the boson transverse momentum. For each sam- 
ple, we fit polynomials to the extracted mean and w^idth of the recoil distributions as functions 
of the boson transverse momentum. The ratios of data to simulation fit-parameters from the 
Z samples are used as scale factors to correct the polynomials parameters of the W simulated 
recoil curves. For each W simulated event, the recoil is replaced with a value drawn from the 
distribution obtained with the corrected parameters corresponding to the W pj. The Sj value 
is calculated by adding back the energy of the W lepton. The energy of the lepton used in 
the calculation is corrected for the energy-scale and resolution effects. Statistical uncertainties 
from the fits are propagated into the Sj distribution as systematic uncertainties. An additional 
systematic uncertainty is included to account for possible differences in the recoil behavior of 
the W and Z bosons. 

The same strategy is followed for the recoil corrections in the electron and muon analyses. As 
an example. Fig. |8] (left) shows the effect of the recoil corrections on the ^j shape for simulated 
events in the electron channel, w^hile Fig. IS] (right) shows the uncertainty from the recoil method 
propagated to the corrected Bj shape of W ^ ev events. The distribution of the residuals, Xf 
is shown at the bottom of each plot, where x is defined as the per-bin difference of the two 
distributions, divided by the corresponding statistical uncertainty. The same definition is used 
throughout this paper. 

The systematic uncertainties on the signal Sj shape are propagated as systematic uncertainties 
on the extracted signal yield through the fitting procedure. Signal shapes are determined for 
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Figure 8: Left: simulated Uj distribution in W — )■ ei/ events before (continuous black line) 
and after (dashed red line) recoil corrections. Right: the uncertainties from the recoil method 
propagated to the corrected Uj shape of W — > ev events (continuous black line, identical to the 
dashed red line on the left-hand side plot) are presented with the red-dashed and blue-dotted 
lines. These two shapes are obtained when the recoil systematic uncertainties are varied by 
one standard deviation. At the bottom of each plot is show^n the distribution of the residuals, 
X, defined as the per-bin difference of the two distributions, divided by the corresponding 
statistical uncertainty. 



the W+ and W separately. 

7.2 Electroweak Backgrounds 

A certain fraction of the events passing the selection criteria for W — > ^i/ are due to other 
EWK processes. Several sources of contamination have been identified. The events with Z — > 
^+£^ (DY background), where one of the two leptons lies beyond the detector acceptance and 
escapes detection, mimic the signature of W — > ^i/ events. Events from Z — > t+t^ and W — > 
TV, with the tau decaying leptonically, have in general a lower-momentum lepton than signal 
events and are strongly suppressed by the minimum pj requirements. 

The Ut shape for the EWK vector boson and tt contributions are evaluated from simulations. 
For the main EWK backgrounds (Z — > C'^l^ and W — > tv), the Uj shape is corrected by means 



of the procedure described in Section 7.1 The ^x shapes are evaluated separately for W+ 

T+v and W^ — ?► t^v. 

A summary of the background fractions in the W — ?► ev and W — ?► //v analyses can be found 
in Table l7| The fractions are similar for the W — ?► ev and W — )■ //v channels, except for the 
DY background which is higher in the W — ?■ ev channel. The difference is mainly due to the 
tighter definition of the DY veto in the W — ?► f/v channel, which is not compensated by the 
larger geometrical acceptance of electrons (I'/l < 2.5) with respect to muons (|?7| < 2.1). 



7.3 Modeling of the QCD Background and W ^> et/ Signal Yield 
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Table 7: Estimated background-to-signal ratios in the W — )■ ev and W — )■ /^v channels. 



Processes 


Bkg. to sig. ratio 


W — )■ ei/ 


W —> }iv 


Z ^ e+e-, ^+^-, T+T- (DY) 

WW+WZ+ZZ 
tt 


7.6% 
3.0% 
0.1% 
0.4% 


4.6% 
3.0% 
0.1% 
0.4% 


Total EWK 


11.2% 


8.1% 



7.3 Modeling of the QCD Background and W — )^ et/ Signal Yield 

Three signal extraction methods are used, which give consistent signal yields. The method 
described in Section |7.3. H is used to extract the final result. 

7.3.1 Modeling the QCD Background Shape with an Analytical Function 

The W — > ei/ signal is extracted using an unbinned maximum likelihood (UML) fit to the Sj 
distribution. 

The shape of the ^j distribution for the QCD background is modeled by a parametric function 
(modified Rayleigh distribution) whose expression is 



/cqd(^t) = ^lexp 



St 



2((7-o + c7i^t)2 



(3) 



The fit to a control sample, defined by inverting the track-cluster matching selection variables 
Afj, A(p, shown in Fig.|9} illustrates the quality of the description of the background shape by the 
parameterized function, including the region of the signal, at high Sj. To study the systematic 
uncertainties associated with the background shape, the resolution term in Eq. dSll was changed 
by introducing an additional QCD shape parameter (72, thus: ctq + (7i^t + (^iSj- 

The free parameters of the UML fit are the QCD background yield, the W signal yield, and the 
background shape parameters ctq and (7i. The following signal yields are obtained: 136328 ± 
386 for the inclusive sample, 81 568 ± 297 for the W+ -^ e+v sample, and 54 760 ± 246 for the 
W^ — > e^v sample. The fit to the inclusive W — > ev sample is displayed in Fig. 10 while the 



fits for the charge-specific channels are displayed in Fig. 11 



The Kolmogorov-Smirnov probabilities for the fits to the charge-specific channels are 0.31 for 
the W+ sample and 0.25 for the W^ sample. Figure 
W sample of the transverse mass, defined as Mj 



12 



shows the distribution for the inclusive 



V2Pt^t(1 
the azimuthal angle between the lepton and the Sj directions. 



cos(A^l^^^)), where A^i^^^ is 



7.3.2 Modeling the QCD Background Shape with a Fixed Distribution 

In this approach the QCD shape is extracted directly from data using a control sample obtained 
by inverting a subset of the requirements used to select the signal. After fixing the shape from 
data, only the normalization is allowed to float in the fit. 

The advantage of this approach is that detector effects, such as anomalous signals in the calo- 
rimeters or dead ECAL towers, are automatically reproduced in the QCD shape, since these 
effects are not affected by the selection inversion used to define the control sample. The track- 
cluster matching variable At] is found to have the smallest correlation with ^j and is therefore 
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Figure 9: Fit to the background-dominated control sample defined by inverting the selection 
on the track-match variables, while maintaining the rest of the signal selection. The blue solid 
line represents the model used to fit the control data sample. This is a Rayleigh function plus a 
floating-yield signal shape that accounts for the signal contamination in the control region. The 
magenta dashed line shows the Rayleigh function alone w^ith its parameters estimated from the 
combined fit. 
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Figure 10: The ^t distribution for the selected W — )■ ev candidates on a linear scale (left) and 
on a logarithmic scale (right). The points w^ith the error bars represent the data. Superimposed 
are the contributions obtained with the fit for QCD background (violet, dark histogram), all 
other backgrounds (orange, medium histogram), and signal plus background (yellowy, light 
histogram). The orange dashed line is the fitted signal contribution. 
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Figure 11: The ^t distributions for the selected W+ (left) and W^ (right) candidates. The 
points with the error bars represent the data. Superimposed are the contributions obtained with 
the fit for QCD background (violet, dark histogram), all other backgrounds (orange, medium 
histogram), and signal plus background (yellowy, light histogram). The orange dashed line is 
the fitted signal contribution. 
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Figure 12: The Mj distribution for the selected W — > ei/ candidates on a linear scale (left) and 
on a logarithmic scale (right). The points with the error bars represent the data. Superimposed 
are the contributions obtained with the fit for QCD background (violet, dark histogram), all 
other backgrounds (orange, medium histogram), and signal plus background (yellowy, light 
histogram). The orange dashed line is the fitted signal contribution. 



chosen as the one to invert in order to remove the signal and obtain the QCD control sample. 
Requirements on isolation and H/E are the same as for the signal selection since these variables 
show significant correlation with I^j. 
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Figure 13: Normalised Sj distribution for QCD and 7+jet simulated events passing the signal 
selection (solid histogram) compared to the normalised distribution for events from all simu- 
lated samples passing the same inverted selection criteria used to obtain the control sample in 
data (dashed histogram). 



The shape of the Sj distribution for QCD and 7+jet simulated events passing the signal selec- 
tion is compared to the Sj distribution for a simulated control sample composed of all simu- 
lated samples (signal and all backgrounds, weighted according to the theoretical production 
cross sections), after applying the same anti-selection as in data (Fig.[l3|. 



The difference in the ^t distributions from the signal and inverted selections is found to be 
predominantly due to two effects, which can be reduced by applying corrections. The first 
effect is due to a large difference in the distribution of the output of a multivariate analysis 
(MVA) used for electron identification in the PF algorithm, between the selected events and 
the control sample. The value of the MVA output determines whether an electron candidate is 
treated by the PF algorithm as a genuine electron, or as a superposition of a charged pion and 
a photon, w^ith track momentum and cluster energy each contributing separately to ^j. The 
control sample contains a higher fraction of electron candidates in the latter category, resulting 
in a bias on the ^j shape. A correction is derived to account for this. The second effect comes 
from the signal contamination in the control sample. The size of the contamination (1.17%) 
is measured from data, using the T&P technique •with Z — )■ e+e^ events, by measuring the 
efficiency for a signal electron to pass the control sample selection. 

The results of the inclusive fit to the Sj distribution w^ith the fixed QCD background shape 



are shown in Fig. 14 the only free parameters in the extended maximum likelihood fit are the 



QCD and signal yields. By applying this second method the following yields are obtained: 
135982 ± 388 (stat.) for the inclusive sample, 81 286 ± 302 (stat.) for the W+ -> e+v sample, 
and 54 703 ± 249 (stat.) for the W^ —?■ e^v sample. The ratios of the inclusive, W+ —> e+v, 
and W^ — ?► e^v yields between this method and the parameterized QCD shape method are 
0.997 ± 0.005, 0.997 ± 0.005, and 0.999 ± 0.005, respectively, considering only the uncorrelated 
systematic uncertainties between the two methods. 

7.3.3 The ABCD Method 

In this method the data are divided into four categories defined by boundaries on ^j and the 
relative tracker isolation, I^^]^/Ej, of the electron candidate. The boundaries of the regions 
are chosen to minimize the overall statistical and systematic uncertainties on the signal yield. 
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Figure 14: Result of the fixed-shape fit to the Sj distribution for all W candidates. The points 
with the error bars represent the data. Superimposed are the results of the maximum likeli- 
hood fit for QCD background (violet, dark histogram), other backgrounds (orange, medium 
histogram), and signal plus background (yellow, light histogram). The orange dashed line (left 
plot) is the fit contribution from signal. 

Values of B'x above and below the boundary of 25 GeV, together with /trk/^T values below the 
boundary of 0.04, define the regions A and B, respectively. Similarly, the regions above and 
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Figure 15: The arrangement of the four categories of events used in the ABCD method. The 
vertical scale indicates increasing values of relative track isolation Ztrk/Ej arid the horizontal 
scale indicates increasing I^j. 



below the Sj boundary for I^^^^/Ej values above 0.04, but below an upper /trk/^x bound of 0.2 
(0.1) for electrons in the EB (EE), define the regions D and C, respectively. There is no upper 
bound for the Sj values. The different regions are shown graphically in Fig. Il5l with region A 
having the greatest signal purity. Combined regions are referred to as 'AB' (for A and B), for 
example. The extracted signal corresponds to the entire ABCD region. 

A system of equations is constructed relating the numbers of observed data events in each of 
the four regions to the numbers of background and signal events, w^ith several parameters to be 
determined from auxiliary measurements or simulations. In this formulation, two parameters, 
/a and /d, relate to the QCD backgrounds and are defined as the ratios of events with a fake 
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electron candidate in the A and D regions to the number in the AB and CD regions, respec- 
tively. The two parameters represent the efficiency with which misidentified electrons pass the 
boundary on ^j dividing AD from BC. If the efficiency for passing the St boundary is largely 
independent of the choice of the boundaries on Ztr^/Ex/ then these two parameters will be ap- 
proximately equal. Assuming /a = /d holds exactly leads to a simplification of the system of 
equations such that all direct dependence of the signal extraction on parameters related to the 
QCD backgrounds is eliminated. For this idealized case there would be no uncertainty on the 
extracted signal yield arising from modeling of QCD backgrounds. Detailed studies of the data 
suggest this assumption holds to a good degree. A residual bias in the extracted signal arising 
from this assumption is estimated directly from the data by studying a control sample obtained 
with inverted quality requirements on the electron candidate, and an appropriate small correc- 
tion to the yield is applied (~0.37%). A systematic uncertainty on the signal yield is derived 
from the uncertainty on this bias correction. This contribution is small and is dominated by the 
uncertainty on signal contamination in the control sample. 

Three other important parameters relate to signal efficiencies: eA arid e^, which are the effi- 
ciencies for signal events in the AB and CD regions, respectively, to pass the Sj boundary, and 
ep, which is the efficiency for the electron candidate of a signal event to pass the boundary on 
relative track isolation dividing the AB region from the CD region under the condition that 
this electron already lies in the ABCD region. The first two of these, ^a arid eo/ are estimated 



from models of the St in signal events using the methods described in Section 7.1 The third 
parameter, ep, is measured from data using the T&P method, described in Section 6.1 and is 
one of the dominant sources of uncertainty on the W boson yield before considering the final 
acceptance corrections. 

Electroweak background contributions are estimated from MC samples with an overall nor- 
malization scaled through an iterative method with the signal yield. 

The extracted yield with respect to the choice of boundaries in relative track isolation and St is 
sensitive to biases in ep and the QCD electron misidentification rate bias correction described 
above, respectively. The yield is very stable with respect to small changes in these selections, 
giving confidence that these important sources of systematic uncertainty are small. 

The foUowng signal yields are obtained: 136 003 it 498 (stat.) for the inclusive sample, 81 525 ± 
385 (stat.) for the W+ -^ e+v sample, and 54356 ± 315 (stat.) for the W" -> e'v sample. 
The ratios of the inclusive, W+ — > e+v, and W^ — > e^v yields between this method and 
the parameterized QCD shape are 0.998 ± 0.007, 0.999 ± 0.007, and 0.993 ± 0.007, respectively, 
considering only the uncorrelated systematic uncertainties betw^een the tw^o methods. 

The results of the three signal extraction methods are summarised in Table IS] 

7.4 Modeling of the QCD Background and W — )^ ;/t/ Signal Yield 

The W — ?> /^t/ analysis is performed using fixed distributions for the St shapes obtained from 
data for the QCD background component and from simulations, after applying proper correc- 
tions, for the signal and the remaining background components. 

Different approaches to signal extraction are considered for W — )■ jiv, as for W — )■ ev. The 
alternative methods do not demonstrate better performance than the use of fixed shapes in the 
W signal fit. Given the lower backgrounds in the muon channel w^ith respect to the electron 
channel, the alternative strategies are not pursued at the same level of detail as in the electron 
case. 
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Table 8: Comparison of W — )■ ev signal extraction methods. The signal yield of each method is 
presented together with its statistical uncertainty. For the fixed shape and the ABCD methods, 
the ratios of the signal yields with the analytical function method are also shown taking into 
account only the uncorrelated systematics between the methods used in the ratios. 



Source 


W ^ ev 


W+ -^ e+v 


W -^ e V 


Analytical fun. 


yield 


136 328 ±386 


81 568 ± 297 


54 760 ±246 


Fixed shape 


yield 
ratio 


135 982 ±388 
0.997 ± 0.005 


81 286 ± 302 
0.997 ± 0.005 


54 703 ±249 
0.999 ± 0.005 


ABCD 


yield 
ratio 


136 003 ±498 
0.998 ± 0.007 


81 525 ± 385 
0.999 ± 0.007 


54356 ±315 
0.993 ± 0.007 



The ?T shape of the QCD background component is obtained from a high-purity QCD sample 
of events that pass the signal selection, except that the isolation requirement is inverted and set 
toJSmb>0-2(Fig.g. 



Simulation studies indicate that this distribution does not accurately reproduce the lij shape 
when muon isolation is required. This is shown in Fig. Il6] (left), where the solid line represents 
the shape for events with an isolated muon and the dashed line the shape obtained by inverting 
the isolation requirement. 
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(right. 



A positive correlation between the isolation variable J^oJ^b ^^'^ ^t is shown in Fig. 
red open circles). This behavior can be parameterized in terms of a linear function ^j oc (1 
a Jjoj^b)' ^^ shown in the same figure. A compensation for the correlation is subsequently made 
by applying a correction of the kind oi Uj = l^j / {\ + oc J^pl^b) *° *^^ events selected by inverting 
the isolation requirement and a new corrected shape is obtained. The agreement of this new 



shape (black points in Fig. 16 left) w^ith the prediction from events with an isolated muon is 
considerably improved. It is also observed that a maximal variation in the correction factor of 
Aa = 0.08 successfully covers the simulation prediction for events with an isolated muon over 



the w^hole ^j interval (shaded area in Fig. 16 left). 

The same positive correlation between ^x and Jjol^b ^^ observed in the data (blue squares in 
Fig. 16 right). A correction Uj = Sj/{1 ± ^I^f^^), with a. ^ 0.2, was applied. The shapes 
obtained in data are shown in Fig. 17 w^here the uncorrected and corrected data shapes from 
events selected by inverting the isolation requirement, together with the simulation expecta- 



tion for events with an isolated muon, are shown. The shaded area in Fig. 17 is bounded by the 
two distributions, obtained using two extreme correction parameters a. ± Aa, with Aa = 0.08, 
as evaluated in simulations. This area is taken as a systematic uncertainty on the QCD back- 
ground shape. 

Several parameterizations for the correction are considered, but the impact on the corrected 
distribution and therefore on the final result is small. Associated uncertainties on the cross 
section and ratios are evaluated as the differences between the fit results obtained with the 
optimal a. value and tw^o extreme cases, a ± Aa. 

The following signal yields are obtained: 140 757 ± 383 for the inclusive sample, 56 666 ± 240 
for the W^ -^ ji^v sample, and 84091 ± 291 for the W+ -^ }i^v sample. 



The Uj distributions are presented in Fig. 18 (full sample) and Fig. n9^ (samples selected by 
the muon charge) superimposed on the individual fitted contributions of the W signal and the 



EWK and QCD backgrounds. Figures 18 and 19 show the Ui distributions for data and fitted 
signal, plus background components. Figure 20 shows the Mj distributions for data and signal. 
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Figure 16: Left: distribution of the corrected Sj for selected events with a non isolated muon 
(black points) superimposed on the distribution of uncorrected Sj for the same events (blue, 
dashed line) and Sj for events with an isolated muon (black, solid histogram). All distributions 
are from simulated QCD events. The shaded area represents the systematic uncertainty due to 
corrections w^ith factors a ± Ace, for Aoc = 0.08. Right: distribution of the average St versus 
^comb ^°^ simulated QCD events (red circles) and for data (blue squares). The high values of St 
in the first two bins in I^f^^ are due to the presence of the W signal events. The superimposed 
lines are linear fits in the range [0.2, 3.0] of Ilf^^,- 
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Figure 17: Distribution of the corrected St for selected events with a nonisolated muon in data 
(black points) superimposed on the uncorrected St distributions for data (blue dashed line) and 
simulated QCD events (black, solid histogram, same as the black, solid histogram in Fig. Il6]|. 
The shaded area represents the systematic uncertainty due to corrections w^ith factors a ± Aa. 
for Aa = 0.08. 
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plus background components, fitted from the Sj spectra. 
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Figure 18: The ^j distribution for the selected W — > /^i/ candidates on a linear scale (left) and 
on a logarithmic scale (right). The points with the error bars represent the data. Superimposed 
are the contributions obtained with the fit for QCD background (violet, dark histogram), all 
other backgrounds (orange, medium histogram), and signal plus background (yellowy, light 
histogram). The black dashed line is the fitted signal contribution. 
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Figure 19: The Sj distributions for the selected W+ (left) and W^ (right) candidates. The 
points with the error bars represent the data. Superimposed are the contributions obtained with 
the fit for QCD background (violet, dark histogram), all other backgrounds (orange, medium 
histogram), and signal plus background (yellow, light histogram). The black dashed line is the 
fitted signal contribution. 



8 The Z -^ i+i- Signal Extraction 

The Z — > £^£^ yield can be obtained by counting the number of selected candidates after 
subtracting the residual background. The Z — )■ i^£^ yield and lepton efficiencies are also de- 
termined using a simultaneous fit to the invariant mass spectra of multiple dilepton categories. 
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Figure 20: The Mj distribution for the selected W — ?► f<v candidates on a Hnear scale (left) and 
on a logarithmic scale (right). The points with the error bars represent the data. Superimposed 
are the contributions obtained with the fit for QCD background (violet, dark histogram), all 
other backgrounds (orange, medium histogram), and signal plus background (yellowy, light 
histogram). The black dashed line is the fitted signal contribution. 

The simultaneous fit deals correctly with correlations in determining the lepton efficiencies 
and the Z yield from the same sample. The Z yield extracted in this way does not need to be 
corrected for efficiency effects in order to determine the cross section, and the statistical un- 
certainty on the Z yield absorbs the uncertainties on the determination of lepton efficiencies 
that would be propagated as systematic uncertainties in the counting analysis. Both methods 
were performed for the Z — > e+e^ analysis, while only the simultaneous fit was used for the 
Z — ?> }i^}i^ analysis after taking into account the results from the previous studies pT|. 

8.1 EWK and QCD Backgrounds 

For the Z — > e+e^ analysis the background contributions from EWK processes Z — > t+t^, tt, 
and diboson production are estimated from the yields of events selected in NLO MC samples 
normalized to the NNLO cross sections and scaled to the considered integrated luminosity. 
They amount to 30.8 ± 0.4 events, where the uncertainty combines the NNLO and luminosity 
uncertainties. Data are used to estimate the background originating from W+jets, 7+jets, and 
QCD multijet events where the selected electrons come from misidentified jets or photons (re- 
ferred to as 'QCD background'). This background contribution is estimated using the distribu- 
tion of the relative track isolation, Jtrk/Ej/ and amounts to 4.9 ± 8.4 (stat.) ± 8.4 (syst.) events. As 
a cross-check, the "same-sign/opposite-sign" method was used, which is based on the signs of 
the charges of the two electron candidates, the measured charge misidentification for electrons 
that pass the nominal selection criteria, and the hypothesis that the QCD background is charge- 
symmetric. The QCD background estimate with this method is 59 ± 17(stat.) ± 160 (syst.) 
events. The two methods are consistent with the presence of negligible QCD background in 
our sample. 

Backgrounds in the Z — )■ f^^f/^ analysis containing two isolated global muons have been es- 
timated with simulations to be very small. This category of dimuon events is defined as the 
"golden" category. The simulation prediction of the smallness of the tt and QCD backgrounds 
was validated with data. First, the selected dimuon sample was enriched with tt events by 
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Table 9: Estimated background-to-signal ratios in the Z — > e+e^ and Z — ?► }i^}i^ (only for 
candidates in the golden category) channels. The QCD background for the Z — ?► e+e^ channel 
has been estimated with data, while all other estimates are based on MC simulations, and their 
corresponding uncertainties are statistical only. 



Processes 


Z -^ e+e sel. 


Z^ }i+ji sel. 


Diboson production 
tt 

Z -^ T+T- 

W+jets 


(0.157 ±0.001)% 
(0.117 ±0.008)% 
(0.080 ± 0.006)% 
(0.010 ± 0.002)% 


(0.158 ±0.001)% 
(0.141 ±0.014)% 
(0.124 ±0.005)% 
(0.008 ±0.002)% 


Total EWK plus tt 


(0.365 ± 0.010)% 


(0.430 ±0.015)% 


QCD 


(0.06 ± 0.14)% 


(0.013 ±0.001)% 


Total background 


(0.42 ± 0.14)% 


(0.444 ±0.015)% 



applying a requirement on Sj, because of the presence of neutrinos in tt events, and an agree- 
ment between data and the simulation prediction was found with the dimuon invariant mass 
requirement inverted, where the residual Z signal is negligible. The QCD component has been 
checked using the same-sign dimuon events and dimuon events with both muons failing the 
isolation requirement, and was found to be in agreement with the simulation predictions. The 
conclusion from the maximum amount of measured data-simulation discrepancy was that the 
uncertainty in the residual background subtraction has a negligible effect on the Z — ?► }i^}i^ 
measured yield. The backgrounds to the Z — )■ fi^ji^ categories having one global and one 
looser muon are significantly larger than in the golden category. Simulation estimates in this 
case are not used for such backgrounds and fits to the dimuon invariant mass distributions are 
performed including parameterized background components, as described in Section [83| 

Backgrounds estimates in the Z — > e+e^ and Z — )■ ji^ji^ analyses are summarized in Tabled] 
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Signal Extraction 



In the following sections the use of a pure Z — )■ e+e^ sample for the determination of the 
residual energy-scale and resolution corrections is first discussed. Then the signal extraction 
with the counting analysis and the simultaneous fit methods are presented. 

8.2.1 Electron Energy Scale 

The lead tungstate crystals of the ECAL are subject to transparency loss during irradiation, 
followed by recovery in periods with no irradiation. The magnitude of the changes to the 
energy response is dependent on instantaneous luminosity and was, at the end of the 2010 
data taking period, up to 1% in the barrel region, and 4% or more in parts of the endcap. The 
changes are monitored continuously by injecting laser light and recording the response. The 
corrections derived from this monitoring are validated by studying the variation of the tt" mass 
peak as a function of time for different regions of the ECAL (using n^ data collected in a special 
calibration stream), and by studying the overall Z — > e+e^ mass peak and w^idth. With the 
current corrections, residual variations of the energy scale w^ith time are at the level of 0.3% in 
the barrel and less than 1% in the endcaps. 

The remaining mean scale correction factors to be applied to the data and the resolution correc- 
tions (smearing) to be applied to the simulated sample are estimated from Z — > e+e^ events. 
Invariant mass distributions for electrons in several tj bins in the EB and EE are derived from 
simulations and compared to data. A simultaneous fit of a Breit-Wigner convolved with a 
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Crystal-Ball function to each Z — > e+e^ n\ass distribution is performed in order to deter- 
mine the energy scale correction factors for the data and the resolution smearing corrections 
for the simulated samples. The energy scale correction factors are below 1% while the resolu- 
tion smearing corrections are below 1% everywhere, with the exception of the transition region 
between the EB and the EE, w^here they reach 2%. Those corrections are propagated in the anal- 
ysis and proper systematic uncertainties for the cross section measurements are estimated as 
discussed in Section 1911 



8.2.2 Counting Analysis 

After energy scale corrections, applied to electron ECAL clusters before any threshold require- 
ment, 10 few^er events (—0.12%) were selected compared to the number of selected events be- 
fore the application of the energy scale corrections. This brings the final Z — > e+e^ sample to 
8442 and, after background subtraction, the Z yield is 8406 ± 92 events. This yield is used for 
the cross section estimation. 

The dielectron invariant mass spectra for the selected sample w^ith the tight selection before 



and after the application of the corrections are show^n in Fig. 21 along with the predicted distri- 
butions. The data and simulation distributions are normalized to account for the difference in 
selection efficiency. 

8.2.3 Simultaneous Fit 

The Z event yield and the electron efficiencies can be extracted from a simultaneous fit. Two 
categories of events are considered: events where both electrons satisfy the tight selection with 
Ej > 25 GeV, and events that consist of one electron with Ej > 25 GeV that passes the tight 
selection, and one ECAL cluster with Ej > 25 GeV that fails the selection, either at the recon- 
struction or electron identification level. 

In each category, a signal-plus-background function is fitted to the observed mass spectrum. 
The signal shape is taken from signal samples simulated with POWHEG at the NLO generator 
level, and is convolved with a Crystal-Ball function modified to include an extra Gaussian on 
the high end tail with floating mean and width. In the first category, the nearly vanishing 
background is fixed to the value reported in Table [9| In the second category of events, the 
background is modeled by an exponential distribution. 

The estimated cross section is 988 ± 10 (stat.) ± 4(syst.) pb. The cross section is in good agree- 
ment with the counting analysis estimate of 992 ±11 (stat.) pb, considering only the statistical 
uncertainty. Both techniques give equivalent results. The counting analysis estimate is used for 
the cross section measurement in the Z — > e^e^ channel. 

8.3 The Z -^ ti~^ti~ Signal Extraction 

The yield of the Z — > ji'^ji^ events is determined from a fit simultaneously w^ith the average 
muon reconstruction efficiencies in the tracker and in the muon detector, the muon trigger effi- 
ciency, as w^ell as the efficiency of the applied isolation requirement. Z — > ^+/^^ candidates are 
obtained as pairs of muon candidates of different types and organized into categories according 
to different requirements: 



• 



7.^^■. a pair of isolated global muons, further split into two samples: 

>-2HI 



• Z^^: each muons associated w^ith an HLT trigger muon; 



• 7}J^: only one of the tw^o muons associated with an HLT trigger muon; 
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Figure 21: Distributions of the dielectron invariar\t mass for the selected Z — ?► e+e^ candi- 
dates on a Hnear scale (top) and on a logarithmic scale (bottom) before (left) and after (right) 
applying energy-scale correction factors. The points with the error bars represent the data. 
Superimposed are the expected distributions from simulations, normalized to an integrated 
luminosity of 36 pb^^. The expected distributions are the Z signal (yellow, light histogram), 
other EWK processes (orange, medium histogram), and tt background (red, dark histogram). 
Backgrounds are negligible and cannot be seen on the linear-scale plots. 
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• Z^s'- one isolated global muon and one isolated stand-alone muon; 

• Zpit: one isolated global muon and one isolated tracker track; 

• Zyj°"^^°: a pair of global muons, of which one is isolated and the other is nonisolated. 

With the exception of the Zl^^'^ category, each global muon must correspond to an HLT trigger 
muon. The five categories are explicitly forced to be mutually exclusive in the event selection: 
if one event falls into the first category it is excluded from the second; if it does not fall into the 
first category and falls into the second, it is excluded from the third, and so on. In this way non- 
overlapping, hence statistically independent, event samples are defined. The expected number 
of events in which more than one dimuon combination is selected is almost negligible. In those 
few cases all possible combinations are considered. 

The five unknowns, the Z yield and four efficiency terms, can be w^ritten in terms of the five 
signal yields in each category as follows: 

t2HLT _ at ^2 ^2 ^2 ^2 



^^^ = ^Z^f(+fi-eHLT^iso%k^sa' (4) 



^fifi - 2Nz^^^+^-eHLT(l -eHLT)eiso^trk^sa' (5) 

N^s = 2Nz^,^+u-eHLTeiso^trk(l -etrk)esa/ (6) 



N^t — 2Nz^^^+^,-eHLTeiso^trk^sa(l — esa)/ (7) 

N]l°i"^^° = 2Nz^^+^-eHLjeiso(l — eiso)%k^sa- (8) 

The various efficiency terms in Eqs. ^ to dsll, the average efficiencies of muon reconstruction in 
the tracker, et^k/ in the muon detector as a stand-alone muon, ega, the average efficiency of the 
isolation requirement, eiso/ and the average trigger efficiency, euu, can be factorized because 
the muon selection factorizes the requirements on the tracker and muon detector quantities 
separately. Neither selection on x^ per degree of freedom nor requirement of the muon recon- 
struction through the tracker-muon algorithm is applied in order to avoid efficiency terms that 
cannot be described as a product of contributions from the tracker and the muon detector. 

The dimuon invariant mass spectra for the five categories are divided into bins of different 
sizes, depending on the number of observed events. The distributions of the dimuon invariant 
mass for the different categories can be written as the sum of a signal peak plus a background 
component. 



Figure 22 shows the dimuon invariant mass spectrum for the Z — )■ /^+f/ golden events on both 



a linear scale and a logarithmic scale, and Figs. 23 and 24 show^ the invariant mass distributions 



for the remaining categories. The spectra are in agreement w^ith the simulation. 

The signal-peak distribution can be considered to be identical in the categories Z^^ and Z^t 
because the momentum resolution in CMS is determined predominantly by the tracker mea- 
surement for muons with pj < 200 GeV. The binned spectrum of the dimuon invariant mass 
in the Z,<^, category, which has the most events of all categories, is taken as shape model for all 
categories but Zpg. The large size of the golden sample ensures that the statistical uncertainty 
of the invariant mass distribution has a negligible effect on the cross section measurement. 
The small presence of background is neglected in this distribution. The uncertainty due to 
this approximation has been evaluated and taken as the systematic uncertainty as described in 
Section |92l 

Because only tracker isolation is used, the shape obtained from golden events can also be used 
to model the Z"""^'*" peak distribution. A requirement on calorimetric isolation would have 
distorted the dimuon invariant mass distribution of events with one nonisolated muon because 
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Figure 22: Distributions of the dimuon invariant mass for the selected Z — > ft^ji^ golden can- 
didates on a linear scale (left) and on a logarithmic scale (right). The points with the error bars 
represent the data. Superimposed are the expected distributions from simulations, normalized 
to an integrated luminosity of 36 pb^^. The expected distributions are the Z signal (yellow, 
light histogram), other EWK processes (orange, medium histogram), and tt background (red, 
dark histogram). Backgrounds are negligible and cannot be seen on the linear-scale plots. 
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Figure 23: Distributions of the dimuon invariant mass for the selected Z^t (left) and Z^s (right) 
candidates. The points with the error bars represent the data. Superimposed are the expected 
distributions from simulations, normalized to an integrated luminosity of 36 pb^^. The ex- 
pected distributions are the Z signal (yellow, light histogram), other EWK processes (orange, 
medium histogram), tt background (red, dark histogram) and QCD background (violet, black 
histogram). 
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Figure 24: Distributions of the dimuon invariant mass for the selected Z^°p^^° candidates. The 
points with the error bars represent the data. Superimposed are the expected distributions from 
simulations, normalized to an integrated luminosity of 36 pb^^. The expected distributions are 
the Z signal (yellow, light histogram), other EWK processes (orange, medium histogram), tt 
background (red, dark histogram), and QCD background (violet, black histogram). 

of FSR, as has been observed both in simulation and data. 

The model of the invariant mass shape for the Z^s category is also derived from golden dimuon 
events. The three-momentum for one of the two muons is taken from only the muon detector 
track fit, in order to emulate a stand-alone muon. To avoid using the same event twice in 
forming the Z^g shape model, the higher-px (lower-px) rnuon is chosen for even (odd) event 
numbers. 

Background shapes are modeled as products of an exponential times a polynomial w^hose de- 
gree depends on the category. Different background models and different binning sizes are 
considered for the categories other than Z^,^ and a systematic uncertainty related to the fitting 
procedure is determined accordingly. 



A simultaneous binned fit based on a Poissonian likelihood 1431 is performed for the different 
categories. Table [lO| reports the signal yield and single-muon efficiencies determined from the 
simultaneous fit and the ratios of the fitted to simulation efficiencies. A goodness-of-fit test 
gives a probability (p-value) of 0.36 for this fit. 

Table 10: Signal yield and efficiencies determined from data with the simultaneous fit, and 
ratios of efficiencies determined from the fit and the simulation. 



Quantity 


Fit results from data 


Data/simulation 


^Z^li+ji- 


13 728 ±121 




eHLT 


0.9203 ± 0.0019 


0.9672 ± 0.0020 


£iso 


0.9813 ± 0.0010 


0.9962 ± 0.0011 


esa 


0.9762 ± 0.0012 


0.9964 ± 0.0013 


etrk 


0.9890 ± 0.0006 


0.9949 ± 0.0007 
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The background in the 2.^^ golden category (of the order of few per mille) was neglected in the 
fit. In order to correct the fitted yield Nz^^+p- for the presence of this background, we subtract 
the small estimated irreducible background fraction. 

A (1.0 ± 0.5)% overall efficiency correction due to the loss of muon events because of trigger 
prefiring is also applied (Section [62] |. 

The estimated cross section is 968 ± 8(stat.) pb. 

9 Systematic Uncertainties 

The largest uncertainty contribution on the measured cross sections is related to the integrated 
luminosity [44J, and amounts to 4%. 



The next most important source of systematic uncertainty is due to the lepton efficiency correc- 
tion factors obtained from the T&P method. In the Z — ?► fi^ji^ analysis, the efficiency uncer- 
tainties are absorbed in the statistical uncertainty of the measurement, via the simultaneous fit 
to the yield and efficiencies. 

Table [TT] shows a summary of systematic uncertainties for the W and Z cross section measure- 



ments. Tables 12 and 13 show a summary of systematic uncertainties for the individual cross 
sections (W+, W^) and the ratios (W+/W^, W/Z). Details of systematic uncertainties for the 
muon and electron channels are described in the following subsections. 

Table 11: Systematic uncertainties in percent for inclusive W and Z cross sections. The "n/a" 
entry means that the source does not apply A common luminosity uncertainty of 4% applies 
to all channels. 



Source 


W ^ et/ 


W ^ }iv 


Z -^ e+e" 


Z -^ ]A+}i- 


Lepton reconstruction & identification 

Trigger prefiring 

Energy/ momentum scale & resolution 

^j scale & resolution 

Background subtraction / modeling 

Trigger changes throughout 2010 


1.4 
n/a 
0.5 
0.3 
0.35 
n/a 


0.9 
0.5 
0.22 
0.2 
0.4 
n/a 


1.8 

n/a 

0.12 

n/a 

0.14 

n/a 


n/a 
0.5 
0.35 
n/a 
0.28 
0.1 


Total experimental 


1.6 


1.1 


1.8 


0.7 


PDF uncertainty for acceptance 
Other theoretical uncertainties 


0.6 
0.7 


0.8 
0.8 


0.9 
1.4 


1.1 
1.6 


Total theoretical 


0.9 


1.1 


1.6 


1.9 


Total (excluding luminosity) 


1.8 


1.6 


2.4 


2.0 



9.1 Electron Channels 

The propagation of statistical and systematic uncertainties on the data /simulation efficiency 
correction factors (p) from the T&P method (reconstruction, identification, and trigger) results 
in uncertainties of 1.4% and 1.8% for the W — ?► ev and Z — ?► e+e^ analyses, respectively. The 
uncertainties on the W+ and W^ cross sections are larger than that for the inclusive W because 
of the larger statistical uncertainty when efficiencies are estimated per charge. The systematic 
uncertainty, which depends on the efficiency under study, is determined by considering alter- 
native signal and background models. The size of the systematic uncertainty is 0.3% for the 
electron selection efficiencies and 1.0% for the electron reconstruction efficiency. The estima- 
tion of the trigger efficiency is considered to be background-free so there is no need to perform 
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Table 12: Systematic uncertainties in percent for individual W cross sections and the ratios in 



nectron channel. A common lummosity i 


mcertamt 


y oi 4 /o a 


jplies to all cross sections. 


Source 


W+(e) 


W- (e) 


W+/W- (e) 


W/Z (e) 


Lepton reconstruction & identification 

Energy scale & resolution 

St scale & resolution 

Background subtraction / modeling 


1.5 
0.5 
0.3 
0.3 


1.5 
0.6 
0.3 
0.5 


1.5 
0.1 
0.1 
0.4 


1.1 
0.2 
0.3 
0.3 


Total experimental 


1.6 


1.7 


1.6 


1.2 


PDF uncertainty for acceptance 
Other theoretical uncertainties 


0.7 
1.0 


1.2 
0.7 


1.6 
1.2 


0.6 
1.2 


Total theoretical 


1.2 


1.4 


2.0 


1.4 


Total (excluding luminosity) 


2.1 


2.2 


2.6 


1.8 



Table 13: Systematic uncertainties in percent for individual W cross sections and ratios in the 
muon channel. A common luminosity uncertainty of 4% applies to all cross sections. 



Source 


W+iji) 


W-{p) 


W+/W- {^) 


W/Z{^) 


Lepton reconstruction & identification 

Trigger prefiring 

Momentum scale & resolution 

St scale & resolution 

Background subtraction / modeling 


0.9 

0.5 

0.19 

0.2 

0.4 


0.9 

0.5 

0.25 

0.2 

0.5 


1.3 



0.06 

0.0 

0.2 


0.9 

0.35 
0.2 
0.4 


Total experimental 


1.1 


1.2 


1.3 


1.1 


PDF uncertainty for acceptance 
Other theoretical uncertainties 


0.9 
0.9 


1.5 
0.8 


1.9 
0.8 


0.9 
1.4 


Total theoretical 


1.3 


1.7 


2.1 


1.6 


Total (excluding luminosity) 


1.7 


2.1 


2.5 


2.0 
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a fit for the signal estimation. Theoretical uncertainties on the corrected efficiencies related to 
the PDF uncertainties and the PDF choice were found to be negligible. 

The electron energy scale has an impact on the Ej distribution for the signal. To study this 
effect, the energy-scale corrections obtained from the shift of the Z mass peak (Section |8.2.1| 
are applied to electrons in the EB and EE in simulation (before the £t requirement) and the 
missing Ej is recomputed. The obtained variations on the signal yield from the UML fit are 
0.5% for the inclusive W, 0.5% for the W+, and 0.6% for the W" samples and 0.1% on the 
W+/W^ ratio. All the charge-related studies (determination of individual W+ and W^ yields 
and W+/W^ ratio and associated systematic uncertainties) include data /simulation charge 
misidentification scale factors, estimated from the fraction of same-sign events in the Z — > e^e^ 
data and simulated samples. 

The energy scale of electrons has an impact on the Z yield because of the Ej > 25 GeV re- 
quirement on the two electrons and the mass w^indow^ requirement. Applying the energy-scale 
corrections mentioned above to the EB and EE electrons and reprocessing the data, the Z yield 
is decreased by 10 events (8452 — > 8442). A systematic uncertainty equal to this decrease of 
0.12% is assigned to the Z signal yield. The energy-scale uncertainty for the W selection is 
included in the systematic uncertainty described in the previous paragraph. There, the system- 
atic uncertainty is larger than that for the Z selection because the energy scale also affects the 
Sj shape used for the signal extraction. The W selection itself is affected by the energy scale at 
the level of 0.12%. 

The Sj shape used in the W fits is also distorted by energy resolution uncertainties; this induces 
a change in the W signal yield by 0.02%. 

The Sj energy scale is affected by our limited knowledge of the intrinsic hadronic recoil re- 
sponse. From the discrepancies found in the data /simulation comparisons (Section |7^, un- 
certainties due to the Sj energy scale are estimated to be 0.3% for inclusive W, W+, and W^ 
yields, and 0.1% for the W+/W" ratio. 

The systematic uncertainties on the background subtraction address the possible difference 
betw^een the true background distribution and the modified Rayleigh function that is used in 
the UML fit. We make the assumption that any such difference can be accounted for by an 



additional dz parameter (defined in Section 7.3 1, which affects the resolution at large values of 



St (below the signal). The value of ax is first determined for three samples: the control sample 
in the data, the control sample in the QCD simulation, and the selected sample in the QCD 
simulation. The values obtained are C72 = 0.0009 GeV"^ 0.0010 GeV"^ and 0.0007 GeV"\ 
respectively for W+ and (T2 = 0.0007 GeV-\ 0.0009 GeV-^ and 0.0008 GeV"^ for W". The 
three values of (72 are then fixed in turn, and ctq and di are set to their values from data to 
generate distributions (of the size of our sample) w^ith the three-parameter function, w^hich 
we then fit with our nominal two-parameter function. The maximal relative difference in the 
yields is quoted as the systematic uncertainty on background subtraction: 0.35% for inclusive 
W, 0.33% for W+, 0.48% for W", and 0.39% for the ratio. 

In the following paragraphs w^e discuss the systematic uncertainties of the fixed shape and the 
ABCD methods w^hich w^ere also explored in order to cross check the extraction of the W — > ei/ 
signal. These uncertainties correspond to the specific methods and are not propagated to the 
final cross section measurement reported in this paper. 

The systematic uncertainties on the background subtraction using fixed-shape distributions are 



summarized in Table 14 The total uncertainty is taken as the sum in quadrature of the values in 



the table, giving 0.40%. The total uncertainty is dominated by the uncertainty of the correction 
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of the signal contamination in the control sample. This uncertainty is evaluated by propagating 
the uncertainty on the measured contamination using the T&P technique. The statistical uncer- 
tainty of this evaluation is calculated using a large number of shapes in which the number of 
events is generated from a Poisson distribution with the mean equal to the number of events 
in the nominal shape. The signal yield under variation of the requirements used to define the 
control sample was also studied and found to be very stable with an RMS spread of 0.12% 
for the range of selections considered. In order to take into account the observation of small 
residual correlations that are not corrected for, an additional systematic uncertainty of 0.35% is 
assigned as a conservative estimate of their size. 

Table 14: Summary of systematic uncertainties for background modeling using fixed-shape 
distributions. 



Source of systematic uncertainty 


Value 


MVA Correction 
Signal contamination 
Statistical fluctuations 
Residual correlations 


0.05% 
0.15% 
0.12% 
0.34% 


Total 


0.40% 



The systematic uncertainties on the signal extraction using the ABCD method are summarized 



in Table 15 The total uncertainty is taken as the sum in quadrature of the individual compo- 
nents listed in the table, and corresponds to 0.7%. The two most important sources of system- 
atic uncertainty arise from the modeling of the signal shape. The largest of these (0.53%) comes 
from the uncertainty on £p, dominated by the statistical uncertainties in the T&P method. Un- 
certainties on the electron energy scale, and hadronic recoil response and resolution affect the 
modeling of the ^j distribution for the signal and together give rise to the second largest un- 
certainty (0.4%). The uncertainty coming from the modeling of electrow^eak backgrounds is 
estimated to be 0.2%. The assumption that the fake electron efficiency to pass the ^j boundary 
is independent of the relative track isolation leads to a small bias for which a correction is ap- 
plied. The uncertainty on this correction gives rise to a very small error on the yield of 0.07%. 
This is dominated by the uncertainty on the signal contamination of the anti-selected sample 
used to estimate the correction. 

Table 15: Summary of systematic uncertainties for the ABCD method. 



Source of systematic uncertainty 


Value 


Signal contamination in bias correction 
EWK backgrounds 
Tag-and-probe 
Sj modeling 


0.07% 
0.20% 
0.53% 
0.40% 


Total 


0.70% 



The QCD background in the Z — ?► e+e^ channel is estimated, as discussed earlier, using the 
shape information of the relative track isolation distribution. The relative uncertainty (approx- 
imately 0.14%) of the total Z yield is used as the systematic uncertainty. 

9.2 Muon Channels 

The total uncertainty of 0.9% (statistical plus systematic) on the correction factors p is used 
as the systematic uncertainty due to muon efficiency (reconstruction, identification, selection. 
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isolation, and trigger) for the W — )■ f/v yield. The systematic uncertainty assigned to the ef- 
ficiencies is evaluated using a large simulated sample including the Z signal and all potential 
backgrounds. Additional uncertainties are evaluated by varying the initial Z preselection cri- 
teria and the mass window to perform the background subtraction fit, and by using alternative 
parameterizations to model the background. The statistical uncertainties on the fit parameters 
describing the background correction are also included. The effect of the uncertainties due to 
the choice of PDFs used in the Z simulation is also studied and found to be negligible. 

The full difference in correction factors for the positively and negatively charged muons, 1.3%, 
is propagated as a systematic uncertainty in the measurement of the W+/W^ cross-section 
ratio. 

A conservative systematic uncertainty of 0.5%, due to the correction for the trigger prefiring in- 
efficiency (Section |6.2| , is assigned to both the Z — > }i^}i^ and W — > /^v cross-section estimates. 

Dedicated studies comparing the peak position and width of the observed Z distribution with 
the expected one indicate a muon momentum scale effect of ~0.25% for 40 GeV muons. In order 
to evaluate the impact on the W cross-section measurement, the fitting procedure with a new 
signal distribution where the muon pj in the simulations is modified according to the observed 
effect, is performed. The difference with respect to the value quoted above is 0.22% for the 
inclusive W sample, 0.19% for W+, and 0.25% for W^, and for the W^"/W~ ratio it reduces to 
0.06%. Muon momentum scale and resolution affect the measurement of the Z — ?► ji^ji^ cross 
section with a 0.35% uncertainty. 

The QCD background shape for the W analysis is tested by applying fits to the Ui spectrum 
with the two extreme Uj shapes, corresponding to the maximal variations of the correction 
factor, oi. The variation in the signal yield with respect to that obtained using the reference 
distribution is 0.4% for the inclusive W sample, 0.4% for W+, 0.5% for W", and 0.2% for the 
W+/W- ratio. 

The recoil modeling in the signal shape is also a potential source of uncertainty. This uncer- 
tainty is estimated by applying the signal shape predicted by the simulation to the fits of the 
^T distribution. The variation in the signal yield with respect to the reference result is 0.2%. 

The systematic uncertainty on the Z — ?► }i^}i^ signal extraction procedure has been evaluated 
as follows. The uncertainty of the fit model is estimated by varying in different ways the back- 
ground models and changing the dimuon mass binning of the various dimuon categories. Half 
of the difference between the maximum and minimum fitted yields across all the tested varia- 
tions is taken as a systematic uncertainty. This amounts to 0.2%. 

The signal shape has been determined assuming that the golden samples are background-free. 
A flat distribution is added as background contribution to the signal shapes and this produces 
a relative change in the fitted Z yield equal to one third of the introduced background fraction. 
An irreducible contamination is known to be present from simulation with the given selection. 
It amounts to less than 0.5%, so a conservative estimate of 0.2% systematic uncertainty due to 
neglecting the background in the signal shapes used for the fit is assigned. Adding those two 
contributions in quadrature, a total systematic uncertainty due to the fit method of 0.28% is 
assigned. 

The stability of the measured Z yields was also checked in the two run periods with different 
trigger thresholds and the corresponding variation of the signal yield of 0.1% is taken as a 
conservative systematic uncertainty. 
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9.3 Theoretical Uncertainties 

The main theoretical uncertainty on the cross section estimation arises from the computation of 
the geometrical and kinematic acceptance of the detector. Uncertainty due to the PDF choice, 
and uncertainties in the PDFs themselves are studied using the full PDF eigenvector set and 
comparing among PDFs provided by the CTEQ, MSTW, and NNPDF groups. For the esti- 
mation of the acceptance uncertainties, we followed the recipe prescribed by the PDF4LHC 
working group ||45| . 

Systematic uncertainties on the acceptances due to the PDF choice are reported in Table [16] 
Here A, denotes the uncertainty (68% confidence level (CL)) within a given set / (i = CTIO [461, 
MSTW08NLO [47J, NNPDF2.1 [48 1). The quantity Agets corresponds to half of the maximum 
difference between the central values of any pair of sets. The final systematic uncertainty (last 
column) considers half of the maximum difference between the extreme values (central values 
plus positive or minus negative uncertainties), again for any pair of the three sets, plus the 
remaining as uncertainties. As can be seen from Table 16 the W^ acceptance uncertainties 
are larger than the W+ ones. This is true for each PDF set as well as for the total assigned 
acceptance uncertainty and reflects the larger d-quark PDF uncertainties with respect to those 
for the u quark. The acceptance estimates obtained using the different PDF sets are summarized 
in Table [iZl 



Table 16: Systematic uncertainties from the PDF choice on estimated acceptances and accep- 
tance correction factors after the analysis selections. 



Quantity 


ACTEQ (%) 


Amstw (%) 


Annpdf (%) 


Asets (%) 


Syst. (%) 


W+ acceptance (e) 


±0.5 


±0.3 


±0.4 


0.2 (NNPDF-MSTW) 


0.7 


W^ acceptance (e) 


±0.9 


±0.5 


±0.7 


0.5 (NNPDF-MSTW) 


1.2 


W acceptance (e) 


±0.5 


±0.3 


±0.4 


0.2 (MSTW-CTEQ) 


0.6 


Z acceptance (e) 


±0.7 


±0.4 


±0.6 


0.3 (NNPDF-MSTW) 


0.9 


W+/W^ correction (e) 


±1.6 


±0.5 


±0.7 


0.7 (NNPDF-MSTW) 


1.6 


W/Z correction (e) 


±0.6 


±0.2 


±0.3 


0.2 (NNPDF-MSTW) 


0.6 


W+ acceptance (ji) 


±0.7 


±0.4 


±0.6 


0.3 (NNPDF-MSTW) 


0.9 


W^ acceptance (ji) 


±1.1 


±0.6 


±0.9 


0.5 (MSTW-CTEQ) 


1.5 


W acceptance (f/) 


±0.7 


±0.4 


±0.6 


0.2 (MSTW-CTEQ) 


0.8 


Z acceptance (f/) 


±1.0 


±0.6 


±0.9 


0.2 (NNPDF-MSTW) 


1.1 


W+/W^ correction (ji) 


±1.9 


±0.6 


±0.9 


0.8 (NNPDF-MSTW) 


1.9 


W/Z correction (pi) 


±0.8 


±0.2 


±0.3 


0.2 (NNPDF-CTEQ) 


0.9 



Table 18 summarizes the different theoretical uncertainties on the acceptance due to ISR and 
NNLO, higher order effects, PDFs, FSR, and missing EWK contributions. 

The baseline MC generator used to simulate the W and Z signals, POWHEG, is accurate up to 
the NLO in perturbative QCD, and up to the leading-logarithmic (LL) order for soft, nonper- 
turbative QCD effects. A description accurate to just beyond the next to next to LL (NNLL) can 
be attained with a resummation procedure Il49ll50t . The ResBos generator [51] implements 
both the resummation and NNLO calculations, which are missing in the baseline generator, 
and its predictions for the W boson pj spectrum show remarkable agreement with pp data at 
^/s = 1.96 TeV 1521. Final state radiation is incorporated in ResBos via PHOTOS ||53l. The 
effect of soft nonperturbative effects, hard higher-order effects, and initial-state radiation (ISR), 
which are not accounted for in the baseline generator, is studied by comparing ResBOS results 
with POWHEG, and the difference is taken as a systematic uncertainty (second column in 



9.3 Theoretical Uncertainties 
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Table 17: Predictions of the central values of the acceptances and the ratios of acceptances for 
various PDF sets. 



Quantity 


CTEQ 


MSTW 


NNPDF 


Aw+(e) 


0.5017 


0.5016 


0.5036 


Aw- (e) 


0.4808 


0.4855 


0.4804 


Aw(e) 


0.4933 


0.4951 


0.4942 


^z(e) 


0.3876 


0.3892 


0.3872 


Aw+(e)/Aw-(e) 


0.9583 


0.9488 


0.9626 


Aw(e)/Az(e) 


0.7857 


0.7853 


0.7880 


Aw+(f) 


0.4594 


0.4587 


0.4617 


Aw-(f) 


0.4471 


0.4519 


0.4472 


Aw(m) 


0.4543 


0.4559 


0.4557 


M^) 


0.3978 


0.3990 


0.3973 


Ay^+{^)/Ay^-{^) 


0.9732 


0.9614 


0.9778 


Aw(m)/Az(h) 


0.8756 


0.8761 


0.8796 



Table 18: Uncertainties on acceptances due to theoretical assumptions. The 
tions are due to ISR plus NNLO effects, factorization and renormalization 
tainties, FSR modeling, and EWK corrections. 



different contribu- 
scales, PDF uncer- 



Quantity 


ISR+NNLO 


liR,}iF Scales 


PDF 


FSR 


EWK 


Total 


W+ acceptance (e) 


0.63% 


0.77% 


0.7% 


0.17% 


0.14% 


1.2% 


W^ acceptance (e) 


0.31% 


0.50% 


1.2% 


0.20% 


0.29% 


1.4% 


W acceptance (e) 


0.53% 


0.34% 


0.6% 


0.13% 


0.14% 


0.9% 


Z acceptance (e) 


0.84% 


0.39% 


0.9% 


0.54% 


0.84% 


1.6% 


W+/W^ correction (e) 


0.32% 


1.14% 


1.6% 


0.26% 


0.25% 


2.0% 


W/Z correction (e) 


0.31% 


0.48% 


0.6% 


0.44% 


1.00% 


1.4% 


W+ acceptance (^) 


0.72% 


0.49% 


0.9% 


0.34% 


0.14% 


1.3% 


W^ acceptance (ji) 


0.50% 


0.37% 


1.5% 


0.16% 


0.39% 


1.7% 


W acceptance (ji) 


0.65% 


0.44% 


0.8% 


0.21% 


0.13% 


1.1% 


Z acceptance {ji) 


1.08% 


0.20% 


1.1% 


0.25% 


1.08% 


1.9% 


W^ /W^ correction {}i) 


0.23% 


0.61% 


1.9% 


0.31% 


0.43% 


2.1% 


W/Z correction (//) 


0.43% 


0.38% 


0.9% 


0.27% 


1.22% 


1.6% 



Table nsl. 



Fixed-order cross section calculations depend on the renormalization (/^r) and factorization 
{jip) scales. Higher-order virtual processes influence the W and Z boson momentum and ra- 
pidity distributions. ResBos fixes Hr and fip to the boson mass, so FEWZ ||54ll55l code is used 
to estimate the effect of scale dependence of NNLO calculations that is quoted as a systematic 
uncertainty. The acceptance is computed by varying up and down the renormalization and 
factorization scales within a factor of two, keeping jIr = j-if. Half of the maximum excursion 
range due to this variation is taken as a systematic uncertainty (third column in Table [18]|. The 



PDF uncertainties from Table 16 are reported in the fourth column of Table 18 and added in 
quadrature to the other contributions to determine the total theoretical uncertainties, shown in 
the last column. 

On top of higher-order QCD corrections, the effect of EWK corrections, not fully implemented 
in our baseline MC samples, is estimated using the HORACE generator Ii56ri59j , which imple- 
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merits both FSR and virtual and nonvirtual corrections. Individual effects are separated and 
the final-state effects are then compared to the PYTHIA results, as PYTHIA is used for FSR 
in the POWHEG event generation. While PYTHIA only partially accounts for NLO EWK 
corrections by generating QED ISR and FSR with a parton shower approximation, HORACE 
also implements one-loop virtual corrections and photon emission from W boson. The differ- 
ence between the two generators is taken as a systematic uncertainty (fifth column in Table p^. 
Moreover, FSR is simulated beyond the single-photon emission in HORACE using the parton 
shower method. The difference due to FSR in HORACE and PYTHIA is taken as a contribution 
to the systematic uncertainty (sixth column in Table [l8]|. 



10 Results 

The results for the electron and muon channels are presented separately. Assuming lepton 
universality, w^e combine our measurements in the different lepton decay modes. The electron 
and muon channels are combined by calculating an average value weighted by the combined 
statistical and systematic uncertainties, taking into account the correlated uncertainties. For 
the cross-section measurements, correlations are only numerically relevant for theoretical un- 
certainties, including the PDF uncertainties on the acceptance values. For the cross section 
ratio measurements, the correlations of lepton efficiencies are taken into account in each lepton 
channel. In the combination of lepton channels, fully correlated theoretical uncertainties are 
assumed for the acceptance factor, with other uncertainties assumed uncorrelated. The lumi- 
nosity uncertainty cancels exactly in the cross-section ratios. 

The NNLO predictions of the total cross sections and their ratios were estimated using FEWZ 
and the MSTW 2008 PDF. The uncertainties, at 68% CL, include contributions from the strong 
coupling DCs Il60l |6T1I , the choice of heavy quark masses (charm and bottom quarks) [62] as 
w^ell as neglected higher-order corrections beyond NNLO, by allowing the renormalization 



and factorization scales to vary in a similar way to that described in Section 9.3 

The following cross sections for inclusive W production are measured: 

a (pp -^ WX) X S (W ^ ev) = 10.48 ± 0.03 (stat.) ± 0.16 (syst.) ± 0.09 (th) ± 0.42 (lumi.) nb, 
a (pp -^ WX) X B{W -^ ^v) = 10.18 ± 0.03 (stat.) ± 0.12 (syst.) ± 0.11 (th) ± 0.41 (lumi.) nb, 
a (pp -^ WX) xB{W -^ £v) = 10.30 ± 0.02 (stat.) ± 0.10 (syst.) ± 0.10 (th) ± 0.41 (lumi.) nb. 

The corresponding NNLO prediction is 10.44 ± 0.27 nb. The results for charge-specific W pro- 
duction are 

a (pp -^ W+X) X B (W+ -^ e+i/) = 6.15 ± 0.02 (stat.) ± 0.10 (syst.) ± 0.07 (th) ± 0.25 (lumi.) nb, 
a (pp -> W+X) X B (W+ ->■ }i+v) = 5.98 ± 0.02 (stat.) ± 0.07 (syst.) ± 0.08 (th) ± 0.24 (lumi.) nb, 
a (pp ->■ W+X) X B (W+ -^ e+v) = 6.04 ± 0.02 (stat.) ± 0.06 (syst.) ± 0.08 (th) ± 0.24 (lumi.) nb, 

and 

a{pp->-'W-X) X B{W- -^e-y) = 4.34 ± 0.02 (stat.) ± 0.07 (syst.) ± 0.06 (th) ± 0.17 (lumi.) nb, 
a (pp -^ W"X) X B (W" ->■ }i-v) = 4.20 ± 0.02 (stat.) ± 0.05 (syst.) ± 0.07 (th) ± 0.17 (lumi.) nb, 
a (pp -^ W^X) X B (W" -^ rv) = 4.26 ± 0.01 (stat.) ± 0.04 (syst.) ± 0.07 (th) ± 0.17 (lumi.) nb. 

The NNLO predictions for these cross sections are 6.15 ± 0.17 nb for W+ and 4.29 ± 0.11 nb 
for W^ . The following cross sections for inclusive Z production are measured: 

a (pp ^ ZX) X B (Z ^ e+e") = 0.992 ± 0.011 (stat.) ± 0.018 (syst.) ± 0.016 (th) ± 0.040 (lumi.) nb, 
a (pp ^ ZX) X S (Z ^ li'^li^) = 0.968 ± 0.008 (stat.) ± 0.007 (syst.) ± 0.018 (th) ± 0.039 (lumi.) nb, 
a (pp -^ ZX) X S (Z ^ l+r) = 0.974 ± 0.007 (stat.) ± 0.007 (syst.) ± 0.018 (th) ± 0.039 (lumi.) nb. 
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The reported Z cross sections correspond to the invariant mass range 60 < nif+f- < 120 GeV, 
and are corrected for the kinematic acceptance but not for 7* exchange. The NNLO prediction 
for Z production is 0.97 ± 0.03 nb. 

The ratio of cross sections for W and Z production is 






Nz ew Aw 



where Az and Aw are the acceptances for Z and W selections, respectively. The two different 
decay channels are combined by assuming fully correlated uncertainties for the acceptance 
factors, with other uncertainties assumed uncorrelated. The resulting ratios are: 



C7(pp -^ WX) X 6(W -^ ev) 



(7(pp- 
a{pp 



ZX) X B{Z - 
WX) X B(W 



e+e-) 



(7(pp -^ ZX) X B{Z -^ ^+^-) 

ajpp -^ WX) X g(W ^ ev) 

(7(pp^ZX) xB{Z^ £+£-) 



= 10.56 ± 0.12 (stat.) ± 0.12 (syst.) ± 0.15 (th.), 
= 10.52 ± 0.09 (stat.) ± 0.10 (syst.) ± 0.17 (th.), 
= 10.54 ± 0.07 (stat.) ± 0.08 (syst.) ± 0.16 (th.). 



The NNLO prediction for this ratio is 10.74 ± 0.04, in good agreement w^ith the measured value. 
The ratio of cross sections for W+ and W^ production is given by 

(7w+ A/w+ £w- Aw- 



(7w- Nw- £w+ A 



w+ 



where Aw+ and Aw- are the acceptances for W+ and W^, respectively. The two different decay 
channels are combined by assuming fully correlated uncertainties for the acceptance factors, 
with other uncertainties assumed uncorrelated. This results in the measurements: 



1.418 ± 0.008 (stat.) ± 0.022 (syst.) ± 0.029 (th.), 
1.423 ± 0.008 (stat.) ± 0.019 (syst.) ± 0.030 (th.), 
1.421 ± 0.006 (stat.) ± 0.014 (syst.) ± 0.029 (th.). 



(7(pp- 


-^ W+X) X B(W+ 


-^ e+i/) 


(7(pp- 
(7(pp- 


-> W-X) X B{W- 
-> W+X) X B(W+ 




(7(pp- 
^(PP- 


-> W-X) X B(W- 
-^ W+X) X e(W+ 


-^e+v) 


cr {pp ^ W-X) X i3(W- 


-^ e-v) 



The NNLO prediction for this ratio is 1.43 ± 0.01, which agrees with the presented measure- 
ment. 



Summaries of the measurements are given in Figs. 25 26 and 27 illustrating the consistency 
of the measurements in the electron and muon channels, as well as confirming the theoretical 
predictions computed at the NNLO in QCD with state-of-the-art PDF sets. For each reported 
measurement, the statistical error is represented in black and the total experimental uncertainty, 
obtained by adding in quadrature the statistical and systematic uncertainties, in dark blue. 
For the cross-section measurements, the luminosity uncertainty is added to the experimental 
uncertainty, and is represented in green. The dark-yellow vertical line represents the theoretical 
prediction, and the light-yellow vertical band is the theoretical uncertainty, interpreted as a 68% 
confidence interval, as described earlier. 
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10 Results 



The ratios of the measurements to the theoretical predictions are listed in Table 19 and dis 



played in Fig. 28 The experimental uncertainty ("exp.") is computed as the sum in quadra- 
ture of the statistical uncertainty and the systematic uncertainties aside from the luminosity 
uncertainty and the theoretical uncertainties associated with the acceptance. The theoretical 
uncertainty ("th.") is computed by adding in quadrature the theoretical uncertainties of the 
acceptance (or the acceptance ratio) and the NNLO prediction, assuming that they are uncor- 
related. 



36pb'at \ls = 7TeV 



36 pb"' at nIs = 7 TeV 



NNLO, FEWZ+MSTW08 prediction 
[with IU1STW08NNLO 68% CL uncertainty] 

10.44 ± 0.27 nb 



W^ev 



W^HV 



W — > Iv (combined) 



NNLO, FEWZ+l/1STW0a prediction, 60-120 GeV 
[with MSTW08NNLO 68% CL uncertainty] 



0.97 ± 0.03 nb 



Z->ee H 

0.992 ± 0.01 1 .,., ± 0.024 .„., + 0.040„.„, nb 
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Figure 25: Summary of the W and Z production cross section times branching ratio measure- 
ments. Measurements in the electron and muon channels, and combined, are compared to the 
theoretical predictions (yellow band) computed at the NNLO in QCD with recent PDF sets. Sta- 
tistical uncertainties are represented as a black error bars, while the red error bars also include 
systematic uncertainties, and the green error bars also include luminosity uncertainties. 



36pb"'at \(s = 7TeV 



36 pb' at xls = 7 TeV 



NNLO, FEWZ+i«STW08 prediction 
[with M5TW08NNLO 68% CL uncertainty] 

6.15 ± 0.17 nb 
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NNLO, FEWZi^iUISTWOS prediction 
[with MSTW08NNLO 68% CL uncertainty] 

4.29 ±0.11 nb 
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Figure 26: Summary of the W^ and W^ production cross section times branching ratio mea- 
surements. Measurements in the electron and muon channels, and combined, are compared 
to the theoretical predictions computed at the NNLO in QCD with recent PDF sets. Statistical 
uncertainties are negligible in this plot; the red error bars represent systematic uncertainties, 
and the green error bars also include luminosity uncertainties. 



Figure 29 show^s the CMS W and Z cross section measurements together with measurements at 
lower center-of-mass energy hadron colliders. The predicted increase of the cross sections with 
center of mass energy is confirmed by our measurements. 

Table |20] reports the cross sections as measured within the fiducial and kinematic acceptance, 
thereby eliminating the PDF uncertainties from the results. In effect, these uncertainties are 
transferred to the theoretical predictions, allowing for a cleaner separation of experimental and 
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36pb'at \ls = 7TeV 



NNLO, FEWZ+MSTW08 prediction 
[with M5TW08NNLO 68% CL uncertainty] 

10.7410.04 



W -> ev, Z ^ ee 

10.560+ 0.120 .,.,,±0.192,,,, 

W -> |j,v, Z ^ |i|j. 

10.52010.090.,., ±0.197.,., 

W -> Iv, Z -> II (combined) 
10.540± 0.070.,., ±0.179..., 
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Figure 27: Summary of the measurements of the ratios of W to Z and W+ to W^ production 
cross sections. Measurements in the electron and muon channels, and combined, are compared 
to the theoretical predictions computed at the NNLO in QCD with recent PDF sets. Statisti- 
cal uncertainties are represented as a black error bars, while the red error bars also include 
systematic uncertainties. Luminosity uncertainties cancel in the ratios. 
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Figure 28: Summary of ratios of the CMS measurements to the theoretical predictions. The 
experimental uncertainties are represented as black error bars, while the red error bars also 
include the combining of theoretical uncertainties on the predictions and measured quantities. 
The yellow band around the vertical yellow^ line at one represent the luminosity uncertainty 
(4%) that affects the cross-section measurements. 



Table 19: Summary of ratios of CMS measurements to the theoretical predictions. 



Quantity 


Ratio (CMS/Theory) 


a X e(W±) 

a X e(W+) 

axB{W-) 

a X b[z) 

axB{W)/axB{Z) 

axB{W+)/(7xB{W-) 


0.986 ± 0.009 (exp) ± 0.028 (th) [±0.029 (tot)] 
0.982 ± 0.010 (exp) ± 0.030 (th) [±0.031 (tot)] 
0.992 ± 0.010 (exp) ± 0.029 (th) [±0.031 (tot)] 
1.002 ± 0.010 (exp) ± 0.032 (th) [±0.034 (tot)] 
0.981 ± 0.010 (exp) ± 0.015 (th) [±0.018 (tot)] 
0.990 ± 0.011 (exp) ± 0.023 (th) [±0.025 (tot)] 
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10 Results 
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Figure 29: Measurements of mclusive W and Z production cross sections times branching ra- 
tios as a function of center-of-mass energy for CMS and experiments at lower-energy colliders. 
The lines are the NNLO theory predictions. 



theoretical uncertainties. For each channel the fiducial and kinematic acceptance is defined as 
the fraction of events with lepton pj greater than 25 GeV (20 GeV for Z — )■ ji^fi^), including 
no final-state QED radiation, and w^ith pseudorapidity in the range \i]\ < 2.5 for electrons and 
\t]\ < 2.1 for muons. 

10.1 Extractionof iB(W-)-€t/) andr(W) 

The precise value of the ratio of the W and Z cross sections obtained from the combination 
of the measurements in the electron and muon final states can be used to determine the SM 
parameters B(W ^ £v) and r(W). 



The ratio of W and Z cross sections can be written as 

(7(pp -^ WX) B(W 



R 



iv) 



(r(pp^ZX) B{Z^i+e- 



In order to estimate the value of B(W — > £v) the predicted ratio of the W and Z production 
cross sections and the measured value of the B{Z — > £^£^) are needed. The NNLO prediction 
of the ratio, based on the MSTW08 PDFs, is c^w/c^z = 3.34 ± 0.08. The current measured value 
for B{Z^ £+£-) is 0.033658±0.000023 HH. Those values lead to an indirect estimation of 



B(W ^ £v) =0.106 ±0.003, 
in agreement with the measured value, ^(W -> £v) = 0.1080 ± 0.0009 

Using the SM value for the leptonic partial width, r(W ^ £v) = 226.6 ± 0.2 MeV 
indirect measurement of the total F ( W) can be obtained through the formula 



, an 



B{W^ £v) 



F(W -^ £v) 

r(w) 
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Table 20: Summary of production cross section measurements in restricted fiducial and kine- 
matic acceptances. The pj and \rj\ requirements restricting the acceptance for electrons and 
muons, and the resulting acceptance values, are also given. The quoted uncertainties on the 
acceptances (evaluated without FSR effect) are due to the PDF uncertainties. 



Channel 


cr X Bin acceptance A (nb) 


A 


W ^ ev 
W+ -^ e+i/ 
W" — > e^v 
Z -^ e+e" 


5.688 ± 0.016 (stat.) ± 0.096 (syst.) ± 0.228 (lumi.) 
3.404 ± 0.012 (stat.) ± 0.067 (syst.) ± 0.136 (lumi.) 
2.284 ± 0.010 (stat.) ± 0.043 (syst.) ± 0.091 (lumi.) 
0.452 ± 0.005 (stat.) ± 0.010 (syst.) ± 0.018 (lumi.) 


0.543 ± 0.003 
0.554 ± 0.004 
0.527 ± 0.006 
0.456 ± 0.004 


pT>25GeV 

|7/|<2.5 


W — > f/v 

w+ -^ ^+v 


4.736 ± 0.012 (stat.) ± 0.067 (syst.) ± 0.189 (lumi.) 
2.815 ± 0.009 (stat.) ± 0.042 (syst.) ± 0.113 (lumi.) 
1.921 ± 0.008 (stat.) ± 0.027 (syst.) ± 0.077 (lumi.) 


0.465 ± 0.004 
0.471 ± 0.004 
0.457 ± 0.007 


Pt >25GeV 

l'/l<2.1 


z ^ ^+^- 


0.396 ± 0.003 (stat.) ± 0.007 (syst.) ± 0.016 (lumi.) 


0.409 ± 0.005 


Pt >20GeV 

1^1 < 2.1 



Based on the above values we obtain 

r(W) =2144 ±62 MeV. 

The SM prediction is 2093 ± 2 MeV 1651 arid the w^orld average of experimental results is 2085 ± 
42 MeV [63 J. The indirect measurement of r(W) is in good agreement with the world average 
and the theoretical prediction, as well as other published measurements. 



11 Summary 

Measurements of the inclusive W and Z production cross sections have been performed using a 
data sample of pp collision events at ^/s = 7 TeV collected with the CMS detector at the LHC in 
2010 and corresponding to an integrated luminosity of 36 pb^^. The inclusive production cross 
sections of W+ and W^ have been measured separately as well as the ratios of the W+/W^ 
and W/Z production cross sections. All measurements are dominated by systematic uncer- 
tainties, the main uncertainty originating from the integrated luminosity (4%), w^hich cancels 
in the ratios. Experimental systematic uncertainties range from 0.7 to 1.8%, and theoretical un- 
certainties range from 0.9 to 2.1%. The measurement of the W/Z cross-section ratio also leads 
to an indirect determination of r(W), which is in agreement w^ith the current world average. 

The results agree with the ATLAS measurement Il20l and with previous CMS results [[2T|. All 
measurements are consistent with the SM NNLO predictions. 
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